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PREFACE 


The simpler observations in X-ray crystallography can 
be treated on the idea that X-rays are reflected from 
layers of atoms in the crystal, and that the beams 
reflected by successive parallel layers interfere with 
one another. This was sufficient, for example, for 
the interpretation of the early ionization-spectrometer 
measurements, and of the Laue photographs of simple 
substances. When more complex crystals were 
studied some means of exhibiting neatly the informa- 
tion given by the large numbers of spots on moving- 
crystal and moving-crystal moving-film photographs 
were necessary, and the device known as the reciprocal 
lattice was found convenient for this purpose. This 
device is indeed almost a necessity for the interpreta- 
tion of such photographs. The reciprocal lattice is 
not generally familiar, and many books on crystal- 
lography stop short at the point where its use becomes 
desirable. This point is roughly the division between 
one- dimensional and three-dimensional problems in 
X-ray optics. Some of the latter problems are treated 
here : the construction and properties of the reciprocal 
lattice, its relation to the commoner types of X-ray 
diffraction photograph, the representation in reciprocal 
space of crystals that are small, distorted or otherwise 
imperfect, and the effects of the imperfections on 
X-ray photographs. The chapter on the effects of 
thermal agitation on X-ray diffraction is included with 
some misgivings, since the writer has no first-hand 
acquaintance with it. One subject that would natur- 
ally form part of X-ray Optics is omitted : the use of 



X-RAY OPTICS 


vi 

Fourier series in the investigation of crystal structures. 
The omission is partly because of the considerable 
space that would be required to deal with it usefully, 
and partly because, in the no-man’s-land of Crystal- 
lography, structure determination is definitely closer to 
Chemistry than to Physics. 

It is hoped that a satisfactory compromise has been 
struck between the difhiseness of a non-mathematical 
treatment and the uninviting symbolism of some of 
the original papers. The necessary mathematics is 
introduced gradually, the earlier chapters requiring 
little more than elementary trigonometry. The mono- 
graphs X-rays and X-ray Crystallography in this series 
have been kept in mind, and it will be found that there 
is practically no overlapping in treatment and very 
little of topic. 

A. J. C. W. 

UNIVERSITY COLLEGE 
CARDIFF 

September igpS 
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CHAPTER I 


INTRODUCTION 

I. Ideally, a crystal * is an arrangement of atoms that 
repeats itself indefinitely in three dimensions. The repeat- 
ing group of atoms occupies a small parallelepiped, called 
the unit cell. The edges of the unit cell are the crystal 
axes, and are generally represented by the letters a, b, c. 
The parallelepiped is not necessarily rectangular, and the 
angles between the pairs of axes are denoted a, y (figure i). 



The complete crystal is made by the regular packing to- 
gether of these units, and the network formed by the edges 
of the unit cells is called the crystal lattice. In most 
crystals the atoms within the unit cell are arranged accord- 
ing to certain geometrical rules of symmetry that govern 
the positions of equivalent atoms. It can be shown that 
there are exactly two hundred and thirty sets of rules 
possible. They are known as space groups, and are of 
fundamental importance in crystal structure analysis. A 
really perfect crystal would be of infinite size, as the 
periodicity necessarily stops at the crystal boundaries, but 

* Fuller descriptions and more precise definitions of the 
technical terms introduced in this paragraph will be found, for 
example, in The Crystalline State^ by W. Lf. Bragg (Bell, I/ondon, 
i933)> or in X^Ray Crystallography, by R. W. James. 
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for most purposes a crystal of more than a thousandth of 
a centimetre in linear dimensions can be considered perfect. 
The main interest of this book is in the diffraction effects 
produced in X-ray photographs of a crystal in which some- 
thing has gone wrong — its size too small, or the lattice- 
network curved, or the space-group rules only approximately 
obeyed. 

2 . It is convenient to begin by recalling some results 
that can be obtained without the use of the reciprocal 
lattice. The first of these is Bragg’s law. In deriving 
this law it is postulated that a crystal consists of parallel 
planes of atoms whose arrangement repeats at regular 



Fig. z . — Derivation of Bragg’s law and estimation of the angular 

range of reflection. 

intervals d (figure 2 ). Such atomic planes are very evident 
when a model of some simple crystal is examined, and may 
be picked out in many different ways. These planes are 
each supposed to reflect a certain small fraction of the 
incident beam of X-rays, so that the angle 6 between the 
incident rays and the planes is equal to the angle betweer 
the reflected rays and the planes, and so that the incidem 
and reflected rays lie in the same plane. For an arbitrar} 
angle 0 the I'cficcted intensity is small, as there is then nc 
special phase relation between the reflections from succes 
sive planes, and the reflected rays largely cancel one an 
other. For certain angles, however, the extra distano 
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travelled by the ray A'B'C' will be an integral number of 
wave-lengths greater than that travelled by the ray ABC^ 
and the reflected rays reinforce one another. It is easy to 
see that, if this condition is fulfilled for the pair of rays 
ABC and A^B'C', it is fulfilled for all other rays, and a 
strong reflected beam is produced. If BD and BE are the 
perpendiculars from B on A'B' and C'B\ the extra distance 
travelled is BB' + B'E, and reinforcement occurs if this 
distance is wA, where n is any integer and X is the wave- 
length of the monochromatic incident X-rays. The lengths 
of T>B' and B'E are each d sin 0, so the beam is strongly 
reflected when 

2 e^sin 0 =wA, . . . (i) 

which is Bragg’s law. The deviation, the angle between 
the directions of the incident and reflected rays, is twice 
the Bragg angle 0. The reflection with w = i is called the 
first-order reflection, that with /z = a is the second order, 
and so on. 

3 . A simple extension of this argument, suggested to the 
writer by Dr A. R. Stokes, leads to an estimate of the 
angle over which the intensity of reflection is appreciable. 
At the angle given by equation (i) all the planes are reflect- 
ing so that their contributions add, but the contribution of 
the last plane is actually mn wave-lengths out of phase with 
that of the first. If it were mn + i or mn — i wave-lengths 
out of phase the intensity of reflection would be zero, since 
the contribution of the top half of the crystal would be 
cancelled by the contribution of the bottom half. As a 
rough estimate the range of 0 over which the intensity of 
reflection is appreciable may be taken as half the difference 
between the angles giving these phase relationships. If 

ze sin O' = (mn + i)A, 
ze sin 0" ^ (mn ~ i)/L, 

£:(sin O' — sin 0") = A, 

O'-vO" . O' -0" . 

sin - = As 


2 S cos 
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or, since d' and 6" differ very little from 6, 

d' — 0" = A/e cos 

where e = md is the thickness of the crystal. If e is small 
(less than about cm.), {O' - 6") is quite appreciable, 

and the crystal reflects over a small range on either side 
of the angle given by equation (i). The effect is analogous 
to the imperfect resolution given by diffraction gratings of 
few lines. Since the actual deviation of the rays is 2.0, 
the range of deviation over which the intensity is appreciable 
is (O' - 6"), or 

^ A/e cos 0. . . . ( 2 ) 

This argument gives only the order of magnitude of /5, 
but the proportionality factor is not far from unity, and 
it is convenient to define the ‘ apparent particle size ’ of 
a crystal by 

e s X/p cos 0 . . , . ( 3 ) 

The intensity of reflection does not, of course, remain 

constant over an angular range and then drop suddenly 
to zero. A more precise definition of ^ in terms of the 
measured variation of intensity of reflection with angle 
is given in chapter II, and its variation with the shape 
and orientation of the crystal particle is discussed in 
chapter IV. 

4 . If the crystal is larger than 10 -^ cm., reflection may 
still take place over a considerable range if the crystal is 
otherwise imperfect. For example, in a distorted crystal 
the interplanar spacing d will not be constant throughout, 
and different parts of the crystal will reflect at different 
angles. The broadening of the lines on powder photo- 
graphs of an aggregate of such distorted crystals may be 
estimated from equation (i). If the maximum strain in 
the crystal is e, the spacing will vary from d{i + e) to 
d(i — d), and 0 will vary over a range numerically equal 
to 2Bd{dd /dd^, or 2 e tan 0. Since the deviation of the rays 
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is 20, the range of angle over which the reflection is appreci- 
able will be twice this, or 

|0/^4etan0. . . . (4) 

In cases of broadening from this cause it is convenient to 
define the ‘ apparent strain ’ by 

Tj . . . ( 5 ) 

Equation 4 suggests that the apparent strain is about four 
times the maximum strain. The relation between actual 
and apparent strain is treated in chapter VIII, but the 
problem is not in as satisfactory a state as that of the 
broadening due to small crystal size. 

5. The other forms of crystal imperfection do not lend 
themselves to elementary treatment, and it seems best to 
postpone discussion of them to later chapters. 



CHAPTER II 


X-RAY DIFFRACTION AND THE RECIPROCAL 

LATTICE 

I. The derivation of the Bragg law given above is simple, 
but it has certain weaknesses. Aesthetically it is displeas- 
ing, because it is a mixture of geometrical and physical 
optics. The interaction of the X-rays with a single plane 
of atoms is treated by geometrical optics, as simple reflec- 
tion, while the interaction with successive planes is treated 
by physical optics, as interference. It is by no means 
obvious why a single plane of atoms ought to reflect X-rays 
as a mirror reflects light ; we shall see shortly that in fact 
it would not. Also, in complex structures, where there 
are atoms in ‘ general ’ positions in the unit cell (that is, 
atoms whose positions are not determined completely by 
symmetry) the successive planes of atoms postulated are no 
longer evident, and the treatment loses its plausibility. We 
must therefore make a new start on the problem of the 
interaction of X-rays and crystals, using physical optics 
throughout, and treating the whole crystal impartially, with 
no bias toward planes and interplanar spaces. This will 
lead naturally to the important concepts of reciprocal space 
and the reciprocal lattice, which form the basis of modern 
X-ray crystallography. 

2. Suppose first that each unit cell contains a single 
scattering centre. This may be imagined as a single atom, 
or as the centre of gravity of the actual distribution of 
electrons in the unit cell. The fact that all the electrons 
are not at the centre of gravity is an additional complication, 
which fortunately can be treated separately as the ‘ structure 
amplitude ’ (James, pp. 49-53, and below, p. 45), and for 
the present it is better to ignore it. In figure i the dots 
represent the scattering centres in the plane containing the 
a and h axes of the crystal, and AB represents the direction 

6 



X-RAY DIFFRACTION AND RECIPROCAL LATTICE 7 

of a parallel beam of monochromatic X-rays falling on it. 
In general, each scattering centre will diifract rays in all 
directions, the amplitude depending on the relation between 
the direction of the incident ray AB and the direction of 
the diffracted ray BC, but there will be some radiation 
diffracted in every direction. Most of the diffracted radia- 
tion is * coherent ' with the incident beam, that is the phase 
of the vibration of the diffracted ray j5C is fixed by the 
phase of the vibration of the incident ray AB. The total 
amplitude of the radiation diffracted in the direction BC 
will therefore be the sum of the contributions of the in- 
dividual scattering centres, proper attention being paid to 



Fig. I. — ^X-rays difTractcd by a plane of scattering centres. 

the phase relations. If it happens that for every scattering 
centre B, B', B",*" the scattered rays BC^ B'C', 
are in phase the radiation diffracted in the direction BC 
will be of great intensity, since all the contributions add. 
On the other hand, if there is no particular relation between 
the phases of the rays BC, B'a, B"C",-- the radiation will 
be weak, as the contributions of the various scattering 
centi'es will mostly cancel each other out. The condition 
that the ray BC and the ray B'C' should add is clearly 
that B'D + B'E should be an integral multiple of A, where 
BD and BE arc perpendicular to A'B' and B'C' respec- 
tively, and A is the wave-length of the incident X-rays. 
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Our object is to find a simple geometrical way of repre- 
senting the fulfilment of this condition for every scattering 
centre of the crystal. This is best done in stages, a line 
of scattering centres being considered first, then a plane, 
and finally the whole crystal. It is convenient, though 
not necessary, to take the line parallel to one of the crystal 
axes, say a, and the plane including two of the crystal axes, 
say a and b. 

In figure 2, then, OX represents a line of scattering 
centres, at a constant distance a apart. The contributions 
of the scattering centres B and B' will add if the path 
difference B'D + B'E = where h is any integer, positive, 
negative, or zero. (The figure is drawn for h positive ; the 



Fig. 2,. — Determination of the path differences between rays 
diffracted by successive scattering centres in a line. 


reader should convince himself that the condition holds for 
other directions of AB and BC, The lines AB^ BC and 
OX are not necessarily coplanar.) A little consideration 
shows that if the contributions of the adjacent scattering 
centres B and B' add, the contribution for any other scatter- 
ing centre B" in the row will add also, the path difference 
for B and jB" being an integral multiple of that for B and B'. 
The condition that B’D + B'E = hX can be represented as 
follows. From a fixed origin O (figure 3) draw a straight 
line of unit length OP parallel to BA, and from P draw 
a straight line of unit length PQ parallel to BC. Let PP' 
and QQ' be perpendiculars from P and Q on OX, where 
OX is parallel to the line of scattering centres. Then the 
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triangle OP'P (figure 3) similar to the triangle P BD 
(figure 2), since their angles are equal, and so 

OP' OP I, 

WD ~ B'B ~ a 

OP' = B'D/a, 

Similarly P'Q' = B'E/a, 

so that OP' + P'Q' == OQ' = {B'D + B'E)/a 

— hX/a 

is the condition that the rays diffracted by all the scattering 
centres B, P', P",-*- add. Interpreted geometrically, 



Fig. 3 . — Geometrical representation of the path differences. 

OQ' =* hX/a is the condition that Q should lie in a plane 
perpendicular to OX and distant hX/a from O. The 
radiation diffracted by the line of scattering centres has 
therefore its greatest possible intensity if Q lies in one of 
a set of parallel planes perpendicular to the line of scatter- 
ing centres and separated X/a from one another. As Q 
moves away from one of these planes the intensity of the 
diffracted radiation will fall off, and will become inappreci- 
able at a distance of the order oi X/t from the plane, where 
t is the length of the line of scattering centres. This is of 
course small compared with X/a^ and for many purposes 
can be neglected. The falling off in intensity is caused by 
the rays diffracted by the various scattering centres becom- 
ing out of phase with one another, and will be considered 
in detail in the next chapter. 
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3. The representation of the condition for reinforcement 
that has just been derived is satisfactory for diffraction of 
X-rays of a single wave-length 1 , but for some uses the 
dependence of the spacing of the planes on A is inconvenient. 
This may, however, be avoided by a simple device. The 
lengths of OP and PQ were taken as unity, as any other 
choice would have seemed artificial as a starting-point, but 
any convenient length may be substituted. If the length 
i/A is chosen for OP and PQ, OP' becomes B'D/Xa^ 
P'Q' becomes B'E/Xa and OQ' = OP' + P'Q' becomes 
(B:D + B'E)/Xa. The condition for reinforcement is 
therefore that OQ' should be equal to h/a^ where h is 
an integer. In geometrical language this is the condition 
that Q should lie in one of a set of parallel planes per- 
pendicular to the < 3 :-axis of the crystal and spaced at distances 
oi 1 /a (figure 4). 



Fig. 4. — Condition that the path differences should be integral 

multiples of the wave-length. 

4- The condition so far derived relates only to a single 
line of scattering centres OX, However, if it is satisfied 
for one line OX, it is satisfied for any parallel line O'X' 
(figure i), since all arguments have involved only the 
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directions of rays, and not the absolute positions of the 
scattering centres. The line OX of figure i could be 
translated bodily to the position O'X' without affecting 
the directions of any of the rays in succeeding figures. 
The case is different with regard to lines of scattering centres 
not parallel to OX. Application of the same argument to 
a line of scattering centres parallel to the i-axis of the 
crystal gives as the condition for the reinforcement of the 
scattered rays that Q should lie on one, say the of a 
set of parallel planes perpendicular to the 6-axis and spaced 
a distance i/h apart. If the two conditions for reinforce- 
ment are to be satisfied simultaneously Q must lie on the 
intersection of these two sets of planes, that is on one of 
a set of lines perpendicular to the plane containing a and 
b and spaced as shown in figure 5. If this condition is 



Fig. 5. — Condition that the path differences between rays dif- 
fracted by all scattering centres in a plane should be integral 
multiples of the wave-length. 


satisfied all centres in the plane shown in figure i scatter 
in phase, and the diffracted radiation is of maximum 
intensity. This behaviour is different from that assumed 
in the derivation of Bragg’s law given in chapter I. A 
little consideration makes it clear that reflection as if by 
a mirror corresponds only to the vertical line through the 
origin in figure 5. 
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If all the scattering centres of the crystal are to diffract 
in phase, a similar condition for the reinforcement of the 
rays scattered by a line of centres in the direction of the 
^-axis must also be fulfilled. As now seems natural, the 
condition is that Q should lie on one, say the of a set 
of parallel planes perpendicular to the c-axis and spaced 
a distance i/c apart. These will intersect the vertical 
lines of figure 5 in a series of points, forming a new lattice 
(figure 6), and all the unit cells of the crystal will scatter 



Fig. 6. — Formation of the reciprocal lattice, expressing the 
condition that the path differences between rays diffracted 
by all scattering centres in a crystal should be integral 
multiples of the wave-length. 


in phase if Q coincides with one of these points. The 
new lattice is characterized by a new set of axes ; that 
having the direction nearest to a is called that nearest 
d is called d*, that nearest c is called and the appropriate 
angles are called oc’^, y*. The crystal axes and the new 

axes coincide in direction only if the crystal ^xes are ortho- 
gonal, and their relative magnitudes are generally different. 
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There are, however, a number of interesting and curious 
■elations between the lattices, which can be summarized 
n the word reciprocal, and the new lattice, representing 
:he conditions that the diffracted intensity should have a 
Tiaximum, is usually called the reciprocal lattice. 

5 . The geometrical conditions for strong diffraction are 
summed up in the integers h, k, /, generally called the 
■ indices of reflection which specify the point of the 
reciprocal lattice with which Q coincides- It is convenient 
to call the direction of the line OQ joining the origin of 
the reciprocal lattice to the point h, k, I the ‘ hkl direction 
The angle between the incident and the diffracted rays 
depends only on the length of OQ, for in the triangle 
GPQ, the angle OPQ is the angle of deviation, and if this 
is called zd, 

S ^ OQ = OP sin 0 + P0 sin (9 
= 2 sin 0/X, 

X = ( 3 /*S) sin Q. 

If J./S is replaced by d/n, this is of the same form as 
Bragg^s law. The distance OQ should therefore be the 
reciprocal of some interplanar spacing in the crystal 
lattice. 

This, or at any rate the converse of it, that for every 
set of planes with spacing d there are points of the reciprocal 
lattice at distances i/d, z/d, '^/d, ••• from the origin, is easy 
to show. In constructing the reciprocal lattice we con- 
sidered three non-parallel lines of scattering centres, which 
for simplicity we chose parallel to the crystal axes. If we 
choose them with other directions, the reciprocal lattice 
resulting must be the same, as it expresses the geometrical 
conditions necessary for the reinforcement of the scattered 
rays, and the final edifice must be the same, whatever the 
scaffolding used in constructing it. Suppose then that we 
take the first two lines of scattering centres lying in one 
of the planes of the set of spacing d, and take the third 
line FG so that one scattering centre belongs to each plane 
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of the set (figure 7). Let the spacing of the centres in 
this line be d\ The condition for reflection from the 

plane as from a mirror is 
that Q should lie on a 
\ \ ^ V^\ \ N ON' through the 

\ V \ y origin of reciprocal space 

\ \ \ perpendicular to the 

\ \ \ \ plane, and the condition 

V \ \ \ ^ for reinforcement of re- 

\ \ ^ flections from successive 

planes is that Q should 
lie on one of a set of 
N * parallel planes perpen- 

dicular to the line of 
scattering centres and at 
a distance x/d' apart. 
Thus Q must be one of 
the points Qq, ••• 

Fig. 7.— Construction to show shovm in the lower part 
that the reciprocal-lattice of figure 7, where ON' 

representation includes £3 parallel to FN, OG' 

to FG, and the spacing 
of the planes is x/d'. The triangles FJH and OKQ^ are 
similar, since their angles are equal, so 


mk 



did':: x/d ' : OQ., 
OQ^ = 


The distances of the points from O are therefore 

x/dy 'z/dy ••*, so that there are points of the reciprocal lattice 
at these distances from the origin. These points correspond 
to the first, second, third, ••• orders of reflection from the 
planes of spacing but the tendency nowadays is to count 
only first-order reflections, and regard the higher orders 
as being due to different sets of planes, with spacings d/z, 
d/ 2 ,y • * *, which happen to be parallel to the set with spacing d. 
Thus for every set of parallel planes in a crystal, there is 
one point of the reciprocal lattice whose distance from the 
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origin is the reciprocal of the spacing of the planes, and 
the line joining the point to the origin of the reciprocal 
lattice (the hkl direction as defined above) is perpendicular 
to the planes. The planes are denoted by the same triplet 
of numbers that describes the point of the reciprocal lattice, 
and it can be shown that this triplet is the same as the 
Miller indices of the plane with the order of reflection 
included (James, pp. 5, 51). 

6. The geometrical construction of paragraph 4 deter- 
mines the reciprocal lattice, but for calculations it is 
convenient to have 
the reciprocal axes 
a*y b*y and 

angles y* 

given explicitly in 
terms of the crystal 
axes ay by c and 
angles a, /?, y. 

The relations be- 
tween the unit 
cells of the lat- 
tices is shown in 
figure 8. Because 
of the method of 
construction of the 
reciprocal lattice 
<2* is perpendicu- 
lar to b and c, b^ to c and <2, and d* to a and b. The 
reciprocal relations follow : a is perpendicular to b^ 
and c*, b to d*^ and a*, and c to <2* and 6*. Also because 
of the method of construction, the projection of a* on a 
is of length i/a, that of b* on b is of length 1/6, and of 
c* on c is of length i/c. The reciprocal relations follow ; 
if the angle between a and is </r, the projection of a* 
on a is a* cos ijjy so a* cos ifj =« i/ay cos ip = i/aa*, and the 
projection of a on a* is a cos ifj = a/aa^ = i/a^. Similarly 
the projection of b on b^ is 1/6^, and of c on c* is i/d^. 



Fig. 8. — Relations between the axes of 
the crystal lattice and those of the 
reciprocal lattice. Right angles are 
marked- 
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From these relations it is easy to find the lengths of the 
reciprocal axes. Let the area of the face of the parallelepiped 
bounded by b and c be Ay that of the face bounded by 
c and a be B, that of the face bounded by a and b be C, 
the volume of the parallelepiped be C/, and let the cor- 
responding functions of the reciprocal axes be denoted by 
the same letters starred. The volume U is equal to the 
area of the face A multiplied by the perpendicular distance 
to the opposite face, and this distance is the projection of 
a on since a* is perpendicular to A. Thus A(i /a*) = 27, 

or == A/U => be sinoL/Uy . . (i) 

and similarly for the other axes. 

If the axes are at right angles each reciprocal axis is 
the reciprocal of the corresponding crystal axis, so that 
27 == abc = i./a^b*d^ = i/U*. The relation 27 = i/27* is 
in fact true for non-orthogonal axes, as the following pro- 
cedure shows. Choose rectangular axes with OJC along a, 
O Y along b^y and OY perpendicular to these and with the 
general direction of c and The direction cosines of the 
crystal and reciprocal axes with respect to these axes are 
readily found to be 

a : ly Oy Oy 

b : cos y, D cosec (cos a — cos P cos y) cosec p, 
c : cos Py o, sin /5, 

a* : Z) cosec a, (cos a cos P — cos y) cosec a cosec /?, 

— T> cosec oc cot Py 

: o, I, o, 

: o, (cos p cos y — cos a) cosec p cosec y, 

D cosec p cosec y, 

where 

D == \/{i —cos^ oc — cos^ P — cos^ y +2 cos oc cos p cos y). ( 3 ) 

Some of these direction cosines can be written down 
immediately, and the rest follow from known angular 
relations between the axes. The volume of a parallelepiped 
is equal to the product of its edges multiplied by the 
determinant of their direction cosines, so 
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U = abcD^ , .... ( 4 ) 

and C/* = cosec oc cosec/? cosec y. (5) 

Replacing by he sin ol/U^ and so on, gives 

= i/t/, . . (6) 

as it was required to show. Use of this result gives 

a* = A /U = sin 

b* = B/U = sin p/bD, I . . (7) 

— C/U — sin y /cD,j 

From the relation reciprocal to equation i 

a == A^/U^ = b^c"^ sin(x*.abcl) 

= a sin a* sin /? sin y/Dy 

since* = Z) cosec /? cosec y, . . . (8) 

or, since cos^ a* = i — sin^ a*, 

cos a* = (cos P cos y — cos a) cosec ^ cosec y, . (9) 

so that 


sin 

sin / 5 * 
sin 

cos a* 
cos / 5 * 
cos y* 


= D cosec P cosec y, 

= JD cosec y cosec oc, 

= D cosec a cosec /?, 

= (cos P cos y — cos a) cosec /? cosec y, 
= (cos y cos a — cos / 5 ) cosec y cosec oc, 
= (cos a cos p — cos y) cosec oc cosec / 5 . 




Kio) 




7. These results make it possible to write down the 
general expression for the Bragg angle in terms of the 
indices hkl. The length S of the line joining the point 
hkl to the origin of the reciprocal lattice may be obtained 
in the following way. Its components along the reciprocal 
axes are /e6*, Ic* ; let it make angles X, fj,y v with these 
axes. The projection of S on any axis is equal to the sum 
of the projections of its components, so 

cos A == ha* ■¥ kb* cos y* + Ic* cos (i i) 

S cos fi = ha* cos y* + kb* -1- Ic* cos oc*, (12) 

S cos V = ha* cos P* + kb* cos a* + Ic*, . (13) 

and S is equal to the sum of the projections of its 
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components on itself, so 

^ ^ ha* cos X + kh* cos ^ + Ic* cos v, . (14) 

Multiplying (ii) by ha*, (12) by hh*, (13) by /c* and adding 
gives 

h^a*"^ + + l^c*^ + 2,klb*c* cos 

+ ‘2,lhc*a* cos 4 - 'zhka*b* cos y* 
== SQia* cos X + kb* cos ja + Ic* cos v) 

= 52, (15) 

by equation 14. Since iS = 2 sin 6/X, the Bragg angle is 
given by 

4 sin^ 6 /X^ = h’^a*^ + k^b*^ + l^c*^ + 2,klb*c* cos a* 

+ 2,lhc*a* cos / 5 * + zhka*b* cosy*. - (ib) 

This equation may be simplified considerably for crystals 
not of the triclinic system. Expressions for each crystal 
system are given in the Internationale Tabellen zur JBestim’- 
mung von Kristallstrukturen, and many books on X-ray 
analysis. Fdf the cubic system with a^b-=c,<x^p^y==* 
90° it reduces to 

^ sin^ d/X^ = (h^ + + l'^)/a^. . (17) 

8 . One reason for the importance of the reciprocal lattice 
is its close relation with X-ray diffraction photographs. 
X-ray photographs resemble the reciprocal lattice much 
more than they resemble the crystal lattice ; one could 
perhaps say that the photograph is merely a more or less 
distorted reciprocal lattice. For special purposes cameras 
have been designed that reproduce a section of the reciprocal 
lattice with little or no distortion (Buerger, 1944). In 
general, however, the distortion is more or less severe, 
and before proceeding with the development of the 
reciprocal lattice it may be well to describe briefly certain 
types of X-ray camera and the photographs they produce. 

9. In the simplest possible arrangement the crystal and 
the film are both stationary and a narrow beam of practically 
parallel X-rays falls on the crystal. If the beam contains 
a wide range of wave-lengths this gives a so-called Laue 
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photograph (James, p. 19), which is not nowadays of 
primary importance. If the beam is monochromatic and 
the crystal perfect there will be in general no strong diffrac- 
tion maxima. The reciprocal lattice and the line OP 
(paragraph 3) are fixed, and strong diffraction would occur 
only if some point of the reciprocal lattice were at a distance 
i/yl from P, or in other words were to lie on a sphere 
(called the * sphere of reflection ’) with centre at P and 
passing through the origin of the reciprocal lattice. For 
an arbitrary orientation of the crystal this is an unlikely 
coincidence. If the crystal is not perfect, however, the 
region of appreciable scattering is not concentrated in a 
point of the reciprocal lattice, but extends to a greater or 
less distance from it, and it is quite likely that the sphere 
of reflection will intersect one or more scattering regions, 
and give rise to more or less diffuse streaks or spots on 
the photograph. Cameras specially designed for investi- 
gating such diffuse regions in reciprocal space are now in 
use, embodying ingenious devices for obtaining a beam as 
nearly monochromatic as possible. Any continuous radia- 
tion would produce Laue spots and make the interpretation 
more difficult. 

10. If the crystal is oscillated or rotated while the ex- 
posure is being made the reciprocal lattice will oscillate or 
rotate through the same angle, many points of the reciprocal 
lattice will pass in turn through the sphere of reflection, 
and the strong scattering occurring as each point intersects 
the sphere will produce corresponding points on the film. 
In the commonest type of moving-crystal camera the crystal 
is mounted in the centre of a cylindrical * film, adjusted 
so that one axis is vertical, and rotated about this axis 
during exposure to a horizontal beam of X-rays. An 
oscillation with uniform angular velocity over a convenient 
range (5°, 10° or 15°) may be substituted for a complete 
rotation. Suppose for definiteness that the c axis is vertical. 

A flat film may be used, but the * layer-lines ’ are then curved 
(James, pp. 24-7). 
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Then the plane of the reciprocal lattice containing and 
b* (the plane = o) is horizontal, and OP lies in it (figure 9). 

The sphere of reflection 
cuts it in a great circle, 
and cuts the parallel planes 
for which Z = ± i, ±2, 
in circles of progressively 
diminishing size. The 
diffracted rays for Z = o all 
lie in a horizontal plane, 
and strike the film along 
the ‘ equator The 
diffracted rays for 
other values of Z lie 
on right circular cones 
whose angle depends on 
the values of Z and c ; the 
semi-angle is easily seen to 
be arc cos {Ik/c). These 
cones intersect the film in 
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Fig. 9. — Formation of layer lines 
on a moving-crystal photo- 
graph. The X-rays emerge 
from the sphere of reflection 
at O, P is the centre of the 
sphere, and four possible 
positions of Q are shown, 
leading to spots with 
Z = 2, I, o and — 2. As the 
crystal moves the reciprocal 
lattice rotates with it, bring- 
ing other points to the surface 
of the sphere of reflection, 
and producing other spots on 
the layer lines. 


circles whose planes are 
parallel to the equatorial 
plane, so that when the 
film is opened out the 
diffraction spots lie on a 
series of parallel ‘ layer 
lines The length of the 
c axis is readily obtained 
from the spacing of these 
lines and the radius of the 
camera. The Z index of a 
spot is equal to the num- 
ber, counted from the 
equator as zero, of the line 


on which it lies. The 


other indices of the spots are generally obtained by a 
method described by Bernal (1926) which involves a 
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certain amount of trial and error, but a direct method 
has been devised (Bunn, Peiser and Turner- Jones 1944 ; 
Turner- Jones and Bunn 1946). 

11. In the Weissenberg camera a cylindrical screen is 
placed between the crystal and the film, with a slot placed 
so as to allow only the spots of one layer line to pass 
through. The film is moved parallel to its axis as the 
crystal is rotated, so that the spots are all over the film, 
instead of being concentrated on a line. It can be shown 
that they lie on two families of curves, from which the 
pair of indices not determined by the layer line can be 
found directly. The method is used chiefly in structure 
determinations, and will not be further described here 
(Weissenberg, 1924). 

12. The methods mentioned so far require a fairly large 
crystal (say not less than o*oi cm. in linear dimensions), 
and if this is not available the powder method can be used. 
In this the specimen, consisting of a large number of small 
crystals of random orientation, is placed in the centre of 
a cylindrical film and bathed in a narrow beam of approxi- 
mately parallel X-rays. During exposure the specimen is 
rotated to make the orientation still more random. By 
this means nearly all directional effects are averaged out, 
the angle of deviation being the only remaining variable. 
The intensity of the X-rays falling on the film with a 
deviation between 20 and 2(0 + if0) will depend on all the 
intensity that, in the representation in reciprocal space, 
lies between spheres of radii = 2 sin 0 /A and 

S + dS = 2 sin 0 /A + ci(2 sin 0 /A). 

If 7 ( 0 ) dO is the intensity diffracted with a Bragg angle 
between 0 and 0 + dO, and 7 (*S) dS is the intensity between 
concentric spheres with radii S and S + dS^ then except 
for certain geometrical factors, 

7 ( 0 ) dO « I{S) dS 

« (2 cos 0 /X)I{S) dO, 

so that 7 ( 0 ) = (2 cos 0 /A) 7 ( 5 '). . . (18) 
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In favourable cases the spheres will not intersect more than 
two regions of high intensity in reciprocal space (say those 
near h, I and near — h, ~ ky — Z), and the photographs 
will consist of a series of lines. Usually the regions round 
hy ky I and round — h, — ky — Z are related by a centre of 
symmetry, and their superposition on the film causes no 
trouble. In unfavourable cases the spheres will intersect 
regions corresponding to several values of hy ky I and the 
interpretation of the film becomes more complicated, but 
powder photographs are of some use as long as discrete 
lines are visible, even when the lines are due to the super- 
position of several regions of high intensity. 

For purposes of structure analysis the important quantity 
connected with a line is its total intensity, that is /sj'7(0) ddy 
the integral being taken over a sufficient length of film to 
include the whole line. Except for certain geometrical 
factors, this is the intensity that would be scattered by 
one unit cell at the angle d multiplied by the number of 
unit cells in the specimen. We shall be more concerned 
with the shape of the line, I{6)y and especially with its 
breadth, or the range of angle over which the reflection 
extends. Various definitions of this have been used, the 
earliest being the half breadth. If Iq is the maximum 
value of I{d)y and are the values of 0 for which 

I{6) = J/o, the half breadth is defined as 

/ 5 ' => 2(02 - dj). . . . (19) 

A definition of breadth better adapted to calculation is the 
integral breadth, defined by 

^ = zI/Iq. . . . ( 30 ) 

(Xhe factor 3 occurs in these expressions because they were 
originally defined in terms of the angle of deviation, not 
the Bragg angle.) If Ig is the maximum value of I(S), 

= IX/Ig cos 61, . . . ( 21 ) 

so that the apparent particle size is 

s = A//5 cos 6 = 7^/7. . 


. (32) 
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13, In deriving the relation betw^een 1 ( 6 ) and I{S) above 
it was tacitly assumed that the variation of intensity on the 
film depended only on the distribution of intensity in 
reciprocal space. This is never true in practice ; the lack 
of parallelism in the incident beam, the finite length and 
diameter of the specimen, and the range of wave-length 
included in even a crystal-reflected X-ray beam would lead 
to a line of finite breadth even if the specimen were ideal. 
The line profile obtained by photometering the film differs 
from that which the crystal imperfection alone would 
produce, and it is necessary to correct it before using it to 
obtain information about the specimen. Several methods 
tiave been devised for doing this, depending on the use of 
;wo lines, one produced by the specimen under investiga- 
:ion, and the other produced under the same experimental 
conditions by a specimen as nearly ideal as possible. From 
:he two line profiles it is possible to derive the line profile 
iue to specimen imperfection alone, but the process is 
:edious and one has generally to be content with the integral 
breadth only. It is found that the breadth of a line pro- 
iuced by a specimen with imperfections is not simply the 
ium of the breadth due to experimental conditions and the 
breadth due to specimen imperfection. The squares of the 
>readths are more nearly additive, but the exact relation 
iepends on the type of imperfection and the experimental 
irrangement. For routine work the correction procedure 
levised by Jones (1938) is perhaps the best. 
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CHAPTER III 


RECIPROCAL SPACE FOR SMALL CRYSTALS 

I. As was seen in the preceding chapter, the radiation 
diffracted by a crystal is at a maximum when the relation 
between the directions of the incident and diffracted rays 
is such that the point Q, found by drawing a line OP of 
length I /A from the origin of the reciprocal lattice parallel 
to the incident ray reversed, and a line PQ from its terminus 
parallel to the diffracted ray, coincides with one of the 
points of the reciprocal lattice. If the coincidence is not 
exact there is still some diffracted radiation, and the next 
step in the investigation is to find the way in which the 
intensity of the diffracted radiation depends on the position 
of Q. Suppose that Q is nearly at the point h, ky I of the 
reciprocal lattice, its co-ordinates in terms of the reciprocal 
axes being h + Uy k ■¥ ‘V, I + w. The successive scattering 
centres in the line OX of figure 2, chapter II, will no longer 
be scattering exactly in phase, and the bigger u is, that 
is the further Q departs from the plane for which 
B'D + B^E = hXy the greater will be the phase difference 
between the rays scattered by successive centres. If 
u were equal to unity, Q would lie on the plane correspond- 
ing to h + 1, and the phase difference would be one wave- 
length, or 2 . 7 t radians. For smaller values of u the phase 
difference between rays scattered by successive centres 
will be proportionately less, and so equal to ZTtu radians. 
If the length of the line of centres is T the number of them 
is T/Uy and the amplitude of the scattered radiation is the 
resultant of AT" s T /a equal vibrations each d = znu radians 
out of phase with the preceding one. 

2. A graphical method of finding the resultant G oi N 
equal vibrations with successive phase differences d is illus- 
trated in figure i. Each vibration is represented by a line 
of unit length, and the phase differences are represented by 

24 
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drawing each line at an angle S to the preceding one. The 
resultant amplitude is given by the distance between the 
free ends of the first and 
last line. In the present 
case TV is large, and the 
polygon made from the 
individual vibrations may 
be considered as the arc of 
a circle, at any rate as long 
as u is small compared with 
unity. The angle subtended 

by the arc at its centre is _ ^ r xr 

1*1 1 • AT Fig. I. — Summation of JV 

I\Oy and Its length is AT, so out-of-phase vibrations. 

the radius of the circle 

is I /<5. The length of the secant joining its ends, which 
is the amplitude of the resultant, is twice the radius 
multiplied by the sine of half the angle subtended at the 



centre, or 


2 . JV(3. 


-sm. 




More accurately, the angle subtended by the arc of the 
circle is TV<3, but its length is A^.sRarcsin {i/zR), where 
R is the radius of the circle. The radius is therefore 
given by 

zNR arc sin {i/zR) = NdR, 
or JR = 1/2 sin {6/z). 

The exact expression for the reflected amplitude is thus 

sin^Ard/2) 

sin {d/z) 


but the expression (i) is a sufficient approximation for all 
applications made here. 

3 . On substituting the values given above for N and <5 
the resultant amplitude becomes 


sin {tiTu/cl) 



3 


TIU 
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The reflected intensity is proportional to the square of 
this, that is to sin^ (jiTu/a) . . 

* ‘ * ^3; 

For u = o this has the value {T/a)^y and as u increases it 
rapidly decreases, reaching zero for u = a/T. It then 
increases a little, but never amounts to more than five per 
cent of its value for u — o. For large ii it falls to zero as 
Figure 2 shows its behaviour as a function of u. 



Fig. 2. — Resultant of N — T/a out-of-phase vibrations as 

a function of u. 


4. We have so far considered the behaviour of a single 
line of scattering centres. In order to investigate the 
behaviour of a complete crystal it is necessary to know its 
shape. In principle it is possible to work out the scatter- 
ing from a crystal of any shape, and it has indeed been done 
for many (Patterson, 1939)* 'T'o illustrate the problem 
two shapes will be considered hei'e, the parallelepiped and 
the sphere. Others can be worked out on the same prin- 
ciples, but the mathematics becomes complicated. The 
general expression is derived in section 6 below. 

Suppose, then, that the crystal has the form of a par- 
allelepiped, whose edges are parallel to the crystal axes, 
and are of lengths -ST, F, Z. The value of T for all lines 
of scattering centres parallel to the <2-axis will be Xy and 
the dependence of the amplitude of reflection on u will be 

sin i^Xu / a) 

7CU 
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for each such line. The dependence of the amplitude on 
V is found by adding Y /h vibrations, each of amplitude 
sin { 7 iXu/d)/nu^ with a phase difference of znv between 
successive vibrations. By the construction of section 2 
the resultant is 

sin {tzXu/o) sin {uzYv/h) 

Till TtV 



Finally, the dependence of the amplitude on 10 is found by 
adding Z/c vibrations, each of amplitude given by expres- 
sion (5), with a phase difference of znw between successive 
vibrations. The dependence of the diffracted amplitude 
G on w, V, w is therefore 

G(u V w) = sm {wcYv/b) sin {tzZw/c) 

7CU TOO TtW 


The diffracted intensity as a function of position in 
reciprocal space, H{u^ v, zo)y is proportional to the square 
of this, that is 


_ sin^ {tcXu/cl) sin^ (TtYv/h) sin^ (tzZw/c) . . 
’ ’ (nu)^ (nv)^ ^ 


The maximum value of H{u, v, zo), for u = v = zv — o, is 
{XYZ)^/{abc)^. Since XYZ is proportional to the volume 

V of the crystal, and abc to the volume TJ of the unit cell, 
the maximum value of H is the square of the total number 
of unit cells in the crystal. For a given distance from the 
point of the reciprocal lattice H has its greatest values 
along the reciprocal axes, that is, the intensity falls off least 
rapidly in directions perpendicular to the faces of the 
crystal. This is found to be true in general ; each plane 
face of a crystal is represented by a ‘ spike ’ of increased 
intensity in reciprocal space, the direction of the spike 
being perpendicular to the face (von Laue, 1936 ; Ewald, 
1938, 1940). If the crystal is plate-like, with, say, X and 

Y much greater than Z, the intensity of reflection as 
a function of position in reciprocal space is practically 
concentrated on the c*-axis, forming a line perpendicular 
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to the plate and of length inversely proportional to its 
thickness. If the crystal is needle-like, with, say, X and Y 
very much less than the intensity of reflection is con- 
centrated into a plate perpendicular to c*, extending a 
little bit further in the directions of and than in the 

intermediate directions. 
These plate-needle and 
needle-plate relations are 
a further illustration of 
the reciprocal nature of 
the relations between the 
two sets of axes. 

5. For a spherical 
crystal the calculation is 
a little more complex, 
and to simplify it as far 
as possible it will be 
assumed that the axes of the crystal are orthogonal. Sup- 
pose that the diameter of the sphere is ID. The length of 

a line of scattering centres at a distance r = — x'^ 

from the <3:-axis is zx (figure 3), so the amplitude diffracted 
by this line is, as a function of w. 



Fig, 3. — Calculation of the in- 
tensity of reflection from a 
spherical crystal. 


sin {zTtxu/d) 

TIU 



The amplitude diffracted by the whole crystal, for 
^ ^ = o, is got by summing this expression for all lines 

of scattering centres. The number of lines of scattering 
centres passing through unit area perpendicular to a is 
equal to a/XJ, so that the sum is approximately 


^ sm. {zTzxu / a) a ' - 

— ZTtr dr 

2)/2 U 


■f 

ro 

= (za/uU) I sin (ztzxu / a)\/ (JD / z)^ 

J D/2 
rD /2 

= {za/uU) I X sin {zjzxu/d) dx. 


- x^ dV{D/zY - 


( 9 ) 
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The integration is easily performed by parts, and gives as 
the amplitude of the diffracted radiation 


G(u) 


7t sin {p:JDu/a) — (tcDu/o) cos {tcDu/o) 
zU (pzu/d)^ 


(10) 


This has been derived as the amplitude diffracted as a 
function of u for v = zo — o. For a spherical crystal, 
however, there is nothing special about the a-axis, and the 
amplitude of the diffracted radiation will depend only on 
the distance 5 — u/ a of Q from the point of the reciprocal 
lattice. This amplitude G(s) will therefore be given by 
equation (10) with u/a replaced by s. The maximum 
value of Gy being the limit of equation (10) for w = o, is 
7 tD^/ 6 Uy which is, as would be expected, just the total 
number of unit cells in the crystal. At the points of the 
reciprocal lattice all unit cells scatter in phase, whatever 
the crystal shape, and the intensity of reflection is propor- 
tional to the square of the number of unit cells. The total 
intensity of a reflection (that is, the integral of the intensity 
over all the region of reciprocal space for which the 
intensity is appreciable) is proportional to the actual 
number of unit cells, so the ‘ integral breadth ’ in reciprocal 
space of a reflection is N/JV^ = i/iV, whatever the crystal 
shape. The shapes of the regions of high intensity, and 
hence the details of the X-ray photographs, are of course 
dependent on the crystal shape. 

6 . Though it is not needed for applications in future 
chapters it is perhaps worth while to write down the general 
expression for the amplitude of reflection as a function of 
position in reciprocal space. For simplicity the treat- 
ment will be limited to orthogonal axes. The rays scat- 
tered by the unit cell that is unit cells in the a direction, 
TZg in the h direction, and in the c direction away from the 
cell at the centre of gravity of the crystal will be advanced 
in phase with respect to those scattered by the cell at the 
centre of gravity by 


(5, /I — Z7c{n^u + + n.^to) 


(II) 
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radians. The amplitude scattered by the whole crystal 
will be the sum of AT = V/U contributions (where V is the 
volume of the crystal and U the volume of a unit cell) of 
equal magnitudes, aqd phases given by equation (ii). 



These cannot be added by the construction given in para- 
graph 2, since depends on three variables, Wi, TZg, n^, 
but the addition can be done trigonometrically by noting 
that the resultant of a number of lines of lengths jPj, 
Fa, •••at angles of d^, da, •*• with the horizontal is repre- 
sented in magnitude G by 

. . . (12) 

and in phase d by 

tan d = B/A, ..... (13) 

where A F-^ cos di + Fg cos da + *** - (14) 

and jB = Fi sin d^ + Fa sin da 4 - •••. 

In the present case the lengths are all unity, and the phases 
are given by equation (ii), so 

A = the sum of cos 2,jt{njU + n^v + n^w) for all values 
of Wi, Wa, W3, 

and B = the sum of sin iZ7t{n-^ + + n^w) for all values 

of T^a, Wg. 

If the crystal is reasonably large m, w are small compared 
with unity, and the summation may be replaced by integra- 
tion without appreciable error. The number of unit cells 
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in the volume element dx dy dz is dx dy dz/U, and x — n^a, 
y = nj), z = n^c, so 

- 4 = cos ^n(ux / a+vy /h+wz / c) dx dy dz^ (15^ 

B = sin zn(ux/a +uy/b +wz/c) dx dy dz, (16) 

H{u, V, w)—A.^ +B^, ...... (17) 

where the integrations extend throughout the volume of 
the crystal. If the shape of the crystal has a centre of 



Fig. 5. — Summation of the contributions from N unit cells. 

symmetry B vanishes. It will be noticed that if u, v, w 
are replaced by - u, — v, — w, A is unaffected and B is 
merely changed in sign. H{u, v, zd) has therefore a centre 
of symmetry at the point of the reciprocal lattice, even if 
the crystal shape is not centrosymmetric. The mathe- 
matically-inclined reader may verify that equations (15)— 
(17) lead to (8) and (lo) for crystals having the shapes of 
a parallelepiped and sphere respectively. 

7. It will be noticed that none of the expressions giving 
the amplitude or intensity as a function of u, v, w depends 
explicitly on h, k, L The surroundings of each point of 
the reciprocal lattice are therefore identical, except perhaps 
for a factor depending on the ‘ structure amplitude ’. This 
should be noted, as it contrasts with the appearance of the 
reciprocal lattice when the crystal is distorted or contains 
* mistakes *. 
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CHAPTER IV 


POWDER PHOTOGRAPHS OF SMALL CRYSTALS 

I. It has been seen that X-ray diffraction is concentrated 
into small regions of reciprocal space in the immediate 
neighbourhood of the points of the reciprocal lattice. The 
diameter of such a region is of the order of i/t, where t is 
a linear dimension of the crystal, and the distance between 
the points of the lattice is of the order of i/a, where a is 
a lattice parameter. The crystals normally used in moving- 
crystal cameras are not smaller than io“^ cms., and a is 
not normally larger than z x cms., so the regions of 

strong diffraction have a magnitude of less than lO"^ of 
their distance apart. The region of strong diffraction is 
therefore inappreciably small, and the finite size of the 
spots on rotating-crystal photographs is not due to diffrac- 
tion in the ordinary use of the word. Geometrical optics 
suffices to explain the size and shape of the spots. With 
powder photographs, however, there is no lower limit to 
the size of crystals that may be used, and if t is of the 
order of io“® cms. the size of the regions of large intensity 
becomes appreciable in comparison with their distance 
apart, and the lines on the powder pattern are broadened. 
It is necessary therefore to determine how the broadening 
depends on the size and shape of the crystals producing 
the powder pattern. 

2. The intensity diffracted with a Bragg angle between 
d and d + dd on 2i powder photograph is proportional to 
the intensity that, in the representation in reciprocal space, 
lies between concentric spheres of radii S = z sin 6 /X and 
S + dS — z sin 6/1 + (2 cos d/X) dd^ and the relation be- 
tween the line-profile I{d) dO and the total intensity included 
between the spheres I{S) dS is (chapter II, paragraph 12) 

7(19) - (2cos( 9/A)7(5). . . (i) 

32 
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The line-profile on the powder photograph can be got 
therefore if the integral of H(Uy v, w) over the space between 
the spheres can be determined. With expressions for H 
as complicated as equation (7) or the square of equation (10) 
of chapter III this is difficult, but it has been done for 
the second one (Patterson, 1939). Normally, however, 
S is very large compared with the regions over which 
H is appreciable, and no error is introduced by replacing 
the spheres K2 by their tangent planes (figure i). 



Fig. I . — Determination of the intensity of reflection as a function 
of 9 from the intensity as a function of position in reciprocal 
space. 

The integral over the region between the planes is difficult 
enough, and is in fact practicable only if the crystals are 
spherical, or for certain reflections from a parallelepiped. 
For spherical crystals it has been done by Patterson (1939), 
and attempts for other crystal shapes have been made by 
replacing the correct expressions for H by approximate 
ones that arc easier to integrate (von Laue, 1936 ; Jones, 
1938 ; etc.). This procedure is, however, of doubtful 
validity, and gives results that are correct only in order of 
magnitude. 

The direct calculation for spherical crystals is somewhat 
too complex to reproduce here, but it gives the same result 
as we shall obtain by a simpler method. Instead of cal- 
culating the amplitude of reflection from the whole crystal 
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as a function of position in reciprocal space, squaring to 
get tke intensity, and integrating to get the intensity with 
S between S and + dS, we shall take the intensity that 
is scattered with S between S and S + dS by a single 
line of scattering centres in the crystal, and sum it for all 
parallel lines. It is obvious that the second process con- 
tains fewer mathematical steps, and it happens in addition 
that they are easier to carry out. 

3. The problem is to calculate the integral breadth of 
a general reflection (that is, one with indices k, 1 ) from 
a crystal of arbitrary shape. The amplitude diifracted by 
a line of scattering centres has been found to be 


G(u) = 


sin (jtTu/a) 

TtU 



where T is the length of the line, a is the repeat distance 
in the direction of the line, and u is the fraction of the 
distance i/a between parallel planes perpendicular to the 
line by which Q misses the nearest point of the reciprocal 
lattice. The intensity scattered with u between u and 
u du is, then 


tt/ \ j sin {^rc'Tu ^ cl\ j . . 

H{u) du = du. . . (3) 

Suppose that the line of scattering centres is parallel to 
the line joining the origin of the reciprocal lattice to the 
point h, ky /, and that Xy y, z are orthogonal axes in the 
crystal with jc in the direction of the line. The total 
intensity scattered by the crystal with u between u and 
u + du will be got by sumrtiing expressions like equation (3) 
for all such lines of scattering centres in the crystal, or to 
a sufficient approximation by multiplying (3) by the 
number of lines per unit area and integrating with respect 
to y and z throughout the crystal. The number of lines 
of scattering centres passing through unit area perpendicular 
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to the lines is a/U, where U is the volume of the unit 
cell, so 


I{u) du = 


a f fsin^ (tiTu/cl) 


dy dz du. 



This quantity J(w) is very closely related to liS) ; in fact, 
the only difference is that u is expressed as a fraction of 
the distance x /a between parallel planes, and S is an actual 
distance in reciprocal space. By replacement of du by 
a dS and u/a by 5 equation (4) becomes 


7(S) dS = C/-1 r dy dz dS, 


J J 


(5) 


5 being the actual distance by which S misses the point 
of the reciprocal lattice.* 

The maximum intensity is given by equation (5) with 


^ = o, so 

I, ^ U- dy dz, . . (6) 

or, since ^dx = T, 

= U-^lTdx dy dz 

^^u-Ardv, . . . (7) 

J V 

where the integration extends throughout the volume of 
the crystal. The total intensity 7 of a reflection is equal 
to the number of unit cells in the crystal, V/U, so that the 
apparent particle size is 

e = I,/I = F-i [ T dV. . . (8) 

J V 


* In this paragraph it has been tacitly assumed that it is possible 
to pick out a line of scattering centres having the direction of the 
line joining the origin of the reciprocal lattice to the point h, k, L 
For a cubic crystal this can be done accurately, and for other 
types of crystal to any desired degree of approximation. If 
equation (4) holds for directions differing slightly from that 
desired, but enclosing it on all sides, considerations of continuity 
require that it hold for the correct direction also. For a derivation 
in which this difficulty does not arise, see Stokes and Wilson 
(1944). 
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The physical significance of (8) is readily seen : the 
apparent particle size is the volume average of the thick- 
ness of the crystal measured perpendicular to the reflecting 
planes. This is a very elegant result, and carries con- 
viction with it, since it seems intuitively probable that the 
broadening of X-ray reflections due to the smallness of the 
diffracting crystal will depend only on the form and 
orientation of the crystal, and not on its particular internal 
construction. The form 

£ = dydz. . . (9) 

obtained from (6) instead of (7), or that given in equation 
(14) below, is frequently more convenient to evaluate. 

The maximum intensity as a function of d is 


Iq = (2 cos 6 /X)Is 

^ {-2. cos e /XU)[ TdV, 


(10) 


J V 


and the integral breadth of the reflection is 

^ = 2///0 

IV 


/ 




(ii) 


cos el TdV 

V 


4. The apparent particle size is easily calculated for 
simple geometrical shapes. For example, for a sphere it 
will be the same for all reflections, and given by (figure 2) 


r2jt rD/2, 

= dy dz = {6/nD^) (D^ - ^r'^)r dr 

J 0 J 0 


dd 


12 rz)^ 


i6_ 


Z)3L8 

~ 3 "^/ .... (i^) 

where D is the diameter of the sphere. It is usual, how- 
ever, to take p ^ ^V SLS the standard ‘ linear dimension ’ 
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of the crystal, and in terms of this 

4V ^ 

= • * • (^3) 

The apparent particle size s is thus just a little less than 
the actual particle size p. 

For crystals other than spherical e will depend on the 
indices of reflection h, k, 1 . For reflection from a cube 




Fig. 2. — Calculation of the apparent particle size for reflections 
from spherical crystals and for the no reflections from 
cubic crystals. 


face it is clearly p, since p is the cube edge, and for planes 
perpendicular to a face diagonal it is (figure 2) 

n 7?/V2 

{^/zp - zy)^ dy dz 
0 

= 2 V 2^/3 

= 0 * 943 ^* ..... (14) 

For this reflection also the apparent particle size is some- 
what less than the actual. 

5. For the practical calculation of the factor connecting 
s and p for a general hkl reflection the form of equation (9) 
is not so convenient as one that can be obtained from it 
by a simple transformation. In figure 3, let A be the 
crystal, and let B be its ‘ ghost shifted a distance t in 
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the hkl direction. Let V(t) be the volume, shaded in the 
hgiire, common to the cr^^stal and its ghost. Then 

V{t) dy djs 

^lR'I'-t)dyd^s, . . 

where T is the length of a line drawn through the Doint 

V, t: in the hkl direction and bounded by the crvstnl 

and the ghost. For t = o. Vit) is the volume o/the 

co'stal, and I ft) vanishes for a value of t dependintr nn 

the size and shape of the crystal. Suppose^ now lihar 

equation (15) is integrated with respect to t fronTt - o 

to a sufficiently large value to make V(t) vanish Th' 
roay be symbolized by ^ vanisti. This 


rm 

J 0 


I:'''"' ■"' ■ iJK ' £‘ ■"] . (.6) 

Clewh simply r, the thickness of the ctystal 

~ed in the „« hhethon thtough the poin. », „ , 

and t dt is 2, so 
J u 


J ^ ^ (0 dss 


= Vs/z, 


* (^ 7 ) 
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y ghost is shifted a distance i in the opposite direction 

t-v* CO JTimon volume is the same as before, so the integral 
^ j^(t^ from a sufficiently large negative value to ^ = o is 

Zlso Ve/z. Thus 


r v(t) dt + dt, 


(18) 


foi" shortness \ V(t) dt, is equal to Ve, and 

e = F-iJF(0 dt. . 

expression may be used to evaluate the relation 
li^e'tw'oen e and p for the hkl reflection from a cube. Sup- 
J30SO 1:tLat the indices are positive and arranged in the order 
Jk. I • This is merely a matter of naming the axes, 
volume common to the crystal and its ‘ ghost ’ is 
thLeo- ^ rectangular parallelepiped of edges p — h' \t\, 

^ k.' \t\^ P ~ \ A> where h', k\ V are the direction cosines 

tli-O reflection, being h, k, I divided by 
*nae limits of integration are t = ±p/h', and the apparent 
P^a^irtiole size is therefore 


= P 


r+p/h' 

{p -h'\t\)(p -k'\t\){p -i'\t\)dt 

J -p/h' 


_ _^r pyh' -{h' +k' +r)p^/2h'^ +{h'k' +ki[^rh')p^/3h'^ 

p^L 

-h'k'i'p^/^'^ 


(/z^ + + r^)( 6 h^ — 2 .hk + kl - 'zlh) 


6 h^ 


(19) 


TTalole: I gives values of the ratio p/& for a number of 
reflections from crystals whose external form is cubic, 
octalxedral or tetrahedral. It is somewhat modified and 
exteixcied from that given by Stokes and Wilson (1942), 
•wlxeire explicit formulae for these shapes and also for the 
rectangular parallelepiped are given. 
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TABLE I 

The ratio of the true to the apparent particle size for various reflections 
from crystals having the form of a tetrahedron^ a cube or an 
octahedron. For a sphere the ratio is constant at 1-0747 


Reflection 

Tetrahedron 

Cube 

Octahedron 

100 

1*3867 

1*0000 

1*1006 

IIO 

0*9806 

1*0607 

1*0376 

III 

1*2009 

1*1547 

1*1438 

210 

1*2403 

1*0733 

1*1075 

3II 

1*1323 

1*1537 

i*io6i 

221 

1*1556 

1*1429 

1*1185 

310 

i' 3 I 56 

1*0672 

1*1138 

31I 

1*2543 

1*1359 

1*121 1 

320 

1*1538 

1*0698 

1*0902 

321 

1*1119 

1*1394 

1*0955 

410 

1*3453 

1*0583 

1*1123 

322, 

1*1772 

1*1556 

1*1304 

41 1 

1*3074 

1*1174 

1*1207 

331 

1*1135 

1*1262 

1*0963 

421 

1*2104 

1*1324 

1*1133 

333 

I *1826 

1-1513 

1*1334 

430 

1*1094 

1*0667 

1*0786 

431 

I *0878 

1*1240 

1*0835 

510 

1*3597 

1*0506 

I*IIOI 


6. The line-profile, that is the intensity of reflection as 
a function of the angle of reflection, is given by equation i : 


= (2 cos 


where I{S) is the integral of the intensity of reflection in. 
reciprocal space over the region bounded by spheres of 
radii S and S + dS. Equation 5 gives a general expression 


for I(S\ 


I r fsin^ {tzTs) 

u] J {Ttsy 


dy dz. 


but, like the corresponding equation (9) for the apparent 
particle size, this expression is not usually the most con- 
venient to evaluate. One alternative form, corresponding 
to equation (7), is 



^sin^ {tzTs) 
(jiTsy 


dV, 


(20) 
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where the integration extends throughout the volume of 
the crystal, but the most useful form is probably 

J(S) = cos (zTcst) dt, . . (21) 

where V(t) is the volume common to the crystal and its 
‘ ghost ’ shifted a distance t in the hkl direction, and the 
limits of integration are the positive and negative values 
of t for which V(t) vanishes. This may be shown by the 
same method as was used in the preceding paragraph. 
Since 


nt) --iKT-t)dyd^ 


(equation 15), 

j* V(t) cos (ZTTSt) dt = iiri: {T — t) cos {2,7tst) dt 

-m: cos (zTzst) d/ — j* ^ cos {zTist) dt \ dy dz 

sin {znst) 


dy dx 


-IK 

-‘II 

-II 


y,sin {ztcsT) _ jySin (zxcsT) 

Z7CS Z7CS 

I — cos (ztcsT) 


f 


Z7ZS 


dy dz 


(pxs)'^ 
sin 2 {tcsT) 

(nsY 


dy dz 


dy dz. 


( 23 ) 


Since 


J V{t) cos (zTisi) dt is equal to J* V{t) cos {zTcst) dt, 
J V{t) cos {zTtst) ^ J'J 


sin^ (tcsT) 


= UI{S), . 


dy dz 


• (23) 


as it was required to show. 

7. Little attention has so far been given to line -profiles, 
since integral breadths are easier to investigate experi- 
mentally. Equation (21) may, however, be evaluated for 
any given value of V(t). Let r be the value of t for which 
4 
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V{t) vanishes. Then, by repeated integration by parts, 

I{S) = cos (27tst) dt 


= 2,U-^ f V(t) cos {ZTtSt) dt 

J 0 

- i^nst) _ {ZTcsi) 


(zjcsy 


(ZTCSy 


+ 




where the primes denote differentiation with respect to t. 
The first term in the square brackets vanishes at both 
limits, since V(r) is zero. For ordinary crystal shapes, 
such as all those mentioned in paragraph 5, V{t) is a cubic 
in t, so that V'" is a constant, is identically zero 

and the final integral vanishes. The line-profile is therefore 
given by 

I(S) — {ns)-^V'{(S) + cos {z7Tsr)V\r) 

-i( 7 ts)-^ sin ( 277 :^T)F"(T)+J( 7 r 5 )-^{i -cos (2:7r5T)}F'"]. (25) 

It is easy to find I{S) for any particular case by substituting 
the appropriate values for V(t) and t in equation (25). 
For example, if the crystal is a cube of edge p, 

v{t) = (p - h't){p - k'txp - I't) 


and T =p/K (paragraph 5), so 

V\t) = - h'{p - k't){p - Vt) - k'(p - rt)(p - h't) 

- lip - h't)(p - k't), 

V'(o) = -PW +k' +/'), 

V'{r) = -PX^' -k'Xh' -l')/h', 

V"(t) = zh!k\p - I't) + zk'l'ip - h't) + zl'h'ip - k't), 
V"{r) = 2p(h'k' - zk'V + I'h'), 

V'" - - Sh'k'V, 
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The line-profile is therefore 


I{S) = 


_ ir p'‘ 


u\_{utsy 


h'+k'+i 


, {h'-k'){W-V) 


h' 


cos (zjcps/h') 


- T~(.h'k' - zk’V + I'h’) sin {zjips/h') 

I jZ/O ) 


2 h'k'l 

2 ,{ 7 CsY 


I — COS {Z7Cps/h') 


_pWh^+k^+l^r 1 f h+k+l (h-k)(h-l) 


2Uh L*An ^ 


A' 


cos Zip 


>p 


I hk — zkl + Ih , , 3 kl, ,.“1 , 


where ijj = {nps'V + l'^)/h. It will be remembered 
that /z, k, I are here the positive values of the indices 
arranged in the order h > k > L 

The asymptotic value of I{S), that is the value about 
wh;ch it oscillates for large values of is given by the 
first term only of equation (25), 


J(5) 


V'jo) 

z{7tsYU^ 


(27) 


so I{S) falls off as s~^. The coefficient of s~'^ is positive, 
since V(£) has its greatest value for ^ = o, and V'(t) must 
be negative for ^ = o. In fact, it is easily seen that — T^'(o) 
is the area of the projection of the crystal on the reflecting 
plane. Equation (27) may therefore be written 

I{S)r^A/zU( 7 tsY, . . . (28) 

where A is this area. An alternative expression for the 
asymptotic behaviour of I{S) may be got from equation (20). 
For ^ large sin {nTs) undergoes rapid alterations in value 
with changes in T (that is, with change in position in the 
crystal), and in the integration over the volume of the 
crystal sin^ (tcTs) may be replaced by its average value. 
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which is one-half. Equation (20) becomes 


7 ( 5 )^ 


I f dV 



which is an interesting contrast to equation (7). The 
value of I{S) for s small depends on the volume average 
of the thickness of the crystal measured perpendicular to 
the reflecting planes ; its value for s large depends on the 
volume average of the reciprocal of the same quantity. 

It is not to be expected that the small oscillations about 
the mean value caused by the trigonometric terms in the 
expressions for the line-profile would ever be observed. 
Even if all the crystals forming the powder specimen had 
the same shape it is unlikely that they would all have the 
same size, and variations in the size would cause the oscil- 
lations to come at different places and so smooth out the 
line-profile. 
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CHAPTER V 


MISTAKES IN LAYER STRUCTURES : GENERAL 

CALCULATION 

MISTAKES IN THREE DIMENSIONS : POWDER 

PHOTOGRAPHS 

I. So far it has been possible to treat each unit cell of 
a crystal as a single scattering centre, and all unit cells as 
identical. If, however, the unit cells are not identical, it 
is necessary to consider the actual distribution of atoms in 
the cell and the combination of scattered rays of different 
amplitudes. The effect of the distribution of atoms can 
be summed up in a quantity called the ‘ structure ampli- 
tude ’. The amplitude of the X-rays scattered by each 
atom can be calculated as a function of the scattering angle 
if the distribution of electrons in the atom is known. It 
is more or less proportional to the atomic number, and is 
greatest at low angles (James, pp. 65-72). It is generally 
denoted by/, with a subscript to distinguish different atoms. 
Tables giving /as a function of sin 6 /A have been prepared 
for most atoms {Internationale Tahellen^ vol. II, pp. 571—4 ; 
Bragg, 1933, pp. 329-34). If / is known for each atom, 
the effect of their different positions within the unit cell 
can be determined. Suppose that there are atoms of scat- 
tering factors /i, /2, ***, /„ at positions in the unit cell 
3^1, % ; -2:2 ; v where y, are given 

as fractions of the edges of the unit cell. For the hkl re- 
flection the rays scattered at (i, i, i) would be advanced 
in phase by 2.Jt{h k + 1 ) radians with respect to those 
scattered at (o, o, o), so the rays scattered at {x, y, z) are 
advanced in phase by 27z{hx + ky + Iss) radians. The re- 
sultant of the rays scattered by the various atoms can be 
constructed graphically as shown in figure i. Each atom is 
represented by a line of length/ at an angle of zuzQix -^ky +/2r) 

45 
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to the horizontal, and the resultant F is represented in 
magnitude by the length of the line joining the beginning 
of to the end of and in phase <5 by the angle that 



O A 

Fig. I. — Graphical construction for the structure amplitude. 

this line makes with the horizontal. Trigonometrically 
expressed, 

. . . (i) 

cos d — A/ |jP|, sin <3 = B/ iF’j, . . ( 3 ) 

where 

A cos + y ’2 cos 2,7t{hx2+ky2+l:^2) 

+ ••• 4-/^, cos Z7t(hXn + ky^ + : . ( 3 ) 

B =f^ sin iZ7z{hx^+ky-^-\-lzj) + sin znQix^^+ky^+l^^ 

+ *** +/n sin ztzQix^ + ky.^ + . . ( 4 ) 

These expressions may be simplified if the positions of some 
or all of the atoms are related by the symmetry elements of 
the space group. For example, in many cases the origin, 
of the unit cell is a centre of symmetry, and an atom at 
the point {x, y, z) implies another at the point (i — x, x — y, 
I — ^). In the equation for B, therefore, for every term 
/sin ‘ZTcifix + Ay + Iz') there is another 

/sin 37r{A(i -x) +A(i -y) +/(i -sr)} /sin 3?r(Ajc+Ay 
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which cancels it, so that B is zero and d is a multiple of tc, '' 
In such cases it is not really necessary to consider the phase 
angle. Since the lines representing the structure ampli- 
tudes are all horizontal, differing merely in direction, they 
may be added algebraically instead of trigonometrically. 
They are taken as positive if <5 = o, negative if (3 = jr. In 
these cases then F = Ay B ^ o. Even if B does not vanish 
for all reflections it is always possible to choose a new origin 
in the unit cell so that it vanishes for the reflection in which 
one is immediately interested. A shift of origin to the 
point, (f, 77, f) would decrease all the angles in figure i by 
Z7t(h^ + k?] + VQ) ; that is, it would merely rotate the figure 
through this angle. Such a rotation would not change the 
magnitude of Fy but the angle <5 could be reduced to zero 
by a proper choice of the new origin. 

The simplifications produced by all the possible com- 
binations of symmetry elements have been tabulated {Inter- 
nationale TahelleUy vol. I ; Lonsdale, 1936). It may 
happen that for some values of h, ky I both A and B are 
zero, and hence that there is no diffracted beam even though 
the conditions for reflection derived in chapter II are satis- 
fied. The corresponding point of the reciprocal lattice is 
then unaccompanied by any region of X-ray diffraction. 
Such reflections of zero intensity are called absences. They 
are of great help in fixing the symmetry of the crystal, but 
care has to be taken not to confuse absences due to the 
symmetry of the crystal with ‘ accidental * absences, which 
occur when 7^^ for a particular reflection happens to be too 
small to be observed. We are, however, concerned not 
with the determination of structures, but with the reciprocal 
space corresponding to a crystal in which F has different 
values in different unit cells. 

2. Mistakes are particularly frequent among the ‘ layer 
structures In these the crystal can be regarded as con- 
sisting of a number of layers, the forces holding each layer 
together being usually much greater than those holding the 
separate layers together. Examples are some silicates such 
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as the micas, cadmium iodide, graphite, and hexagonal 
cobalt. If the forces between the layers are weak the layers 
may be displaced parallel to themselves, the arrangement 
of atoms within the layer being maintained, but the proper 
sequence of atoms in passing from one layer to the next 
being impaired. Even if the forces between the layers are 
strong, there may be several positions in which one layer 
can fit on to that below it without much difference in energy, 
and so occasional mistakes are possible. Cobalt is an 
example of this. It consists of close-packed planes of 
atoms, such as those shown in figure z. If a plane of atoms 
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Fig. 2 ,. — Possible positions of planes of atoncis in close packing. 


has the position A, the succeeding one can take either the 
position B or the position C. Suppose it takes the posi- 
tion B. The third plane could take either the position C 
or A. If the cobalt is in the hexagonal modification it 
should take the position A, and successive planes should 
maintain the sequence ABABABABA---. As far as the 
atoms in contact are concerned, however, the combination 
AB can be followed equally well by A or C, and it would 
not be surprising if an occasional plane forgot to consult 
the atoms with which it was not in contact, and settled 
down in the wrong place. Photographs of hexagonal cobalt 
do show peculiarities most easily explicable by mistakes of 
this kind. The layer silicates have much more complicated 
structures, but many of them alsa have three (or perhaps 



MISTAKES IN LAYER STRUCTURES 49 

six) possible positions of equilibrium. With some of them 
there is another possibility of mistakes. Layers of water 
or other substances can be adsorbed in the spaces between 
the layers, thus increasing the interlayer spacing. If the 
number of layers of water molecules is the same in every 
interlayer space the crystal remains perfect, but if the num- 
ber of layers of water molecules varies there will be erratic 
variations in the interlayer spacing. 

3. Let us first consider a hypothetical case in which the 
a and h axes lie in the plane of the layer and c is perpen- 
dicular to it, and the mistake consists of a shift of the layer 
by half the h axis. If F is the structure factor of a cell in 
an unshifted layer, the structure factor of a cell in a shifted 
layer is easily seen to be (— )^jP, since the y- coordinate of 
every atom is increased by 1/2, and both A. and B are 
changed in sign for k odd, and are unchanged for k even. 
If we try to calculate the amplitude of reflection by the 
method of chapter III, paragraph 4, we find that the sum- 
mations for lines parallel to the a and h axes in the 
layer go through successfully, leading to the expression 

^ sin {tiXu/cl) sin {jiYv/h) 

r ^ ^ . > 

7ZU TtV 

where is F or {—Y'F depending on whether or not the 
layer is a shifted layer, but the third step, of adding 
the amplitudes diffracted by successive layers to give the 
dependence on w, cannot be taken, since the layers are not 
identical. As we are considering the effect of mistakes on 
reciprocal space, not the effect of small particle size, we 
may suppose that AT, Y and are large, and the result 
therefore means that the reflection is practically concen- 
trated into a region with very little extension in the direc- 
tions of <2* and 6*, but that there may be a considerable 
spread in the direction of We have now to find some 
way of estimating this- 

Suppose that the crystal starts with a layer in the ‘ right * 
position, and that each added layer tends to be in the right 
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position, but that there is a chance a of its going into the 
‘ wrong * position. If it goes into the wrong position, 
every layer after it will go into the wrong position, until 
there is another mistake, which restores the layers to the 
right position. Suppose that is the probability of 
the added layer being in the right place, and the 
probability of its being in the wrong place. Clearly 

Rn + I • . • ( 5 ) 

The probability of the (n + added layer being right 
can be calculated as follows. If the added layer is 
right the (n + will be right unless there is a mistake, 
and if the added layer is wrong the (n + will be 
right only if there is a mistake. The probability of the 
(n + being right is therefore the probability of the 
being right multiplied by the probability of no mistake, 
plus the probability of the being wrong multiplied by 
the probability of there being a mistake. In symbols this 
becomes 

Rn+l ~ 

= (i - a)R^ + a(i - R^) 

= a + (i - 2oc)Rn . . • (6) 

This equation may be solved by trying a solution of the 
form 

Rn=P + Qp^ . . . ( 7 ) 

where P, Q and p are constants to be determined. Sub- 
stitution gives P = -J, p — I — 2a, so 

= i + e(i - . . (8) 

Since Pq is i, the starting layer being correct, Q must 
have the value and 

P71 ^ “ 2a)^} . . (9) 

In the next paragraph the average value of the product of 
the P*s of two layers separated by n interlayer spaces will 
be needed. This is clearly P^ multiplied by the proba- 
bility of the two layers being the same (that is, both right 
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or both wrong), plus multiplied by the probability 

of the two layers being different. A little consideration 
shows that these probabilities are and so that the 
required average is 

F^R^+ (-fF^{i - R^) 

= F^, k even, .... (10) 

= (i — zolYF^, k odd. . . (ii) 

4. Let us now try to find the intensity of reflection for 
w — o, that is at the point where the intensity is concen- 
trated for a perfect crystal. The amplitude is given by 

G = + ^3 + ••* + . . (12) 

where N ^ Z/c and F^^ the structure factor of the 
layer, is equal to F for k even, and is either F or — F for 
k odd. If k is even therefore G = NF = -ZF/c, just as for 
a perfect crystal. If k is odd the terms in equation (12) 
will be positive for a while, then negative for a while, 
positive again, and so on, so that the sum is indetermihate. 
It is, however, possible to find its probable value, or rather 
its probable value without regard to sign. Depending on 
how the mistakes occur, G could have any value between 
— NF and + NF, and its average value will be zero, but 
usually it will have some finite sum either positive or nega- 
tive. In calculating the X-ray intensity the sign does not 
matter, and it is the average value of that is wanted. 
Squaring equation (12) gives 

G2 = i?2 + iv’l + ... -i- 

+ 2(Fi>2 + F^F^ + F3F4 + 4. 

+ 2 .{F -^F 3 4 - F qF 4 4 * F 3F 5 4 - • * * + F 2F . 

_j_ ... 

+ ^F^F^, . . . . . * (^3) 

where the terms have been arranged so that in the first 
line the two F’s multiplied together refer to the same layer, 
in the second they belong to adjacent layers, in the third 
there* is one intervening layer, and so on. Suppose that 
is the average value of the product of the F’s of two 
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layers separated by m interlayer spaces. Then 

= NJo + 2(iV - i)/i + 2 (-^ - 2)/2 + — + zJiir-i . (r4) 

For the particular mistake considered in the last paragraph 
was found to be 

= (r - . . . (Ii) 

so 

Q2 ^ F^\2,{N {N - i)(i — 2a) + {N - 2)(i - 2 <x,y 

+ ••• + (i - 2a)^'-i } ” iV]. . (15) 

The expression within the curly brackets is a combined 
arithmetic-geometric series, and it is not difficult to find 
its sum. In the present case, however, N is large and 
(i — 2a)”^ decreases rapidly, so it may for practical pur- 
poses be regarded as geometric series of the form N{i — 2a)^ 
summed to infinity, and 

= 2NF^{i + (i — 3a) + (i - 2a)2 + •••} 

= zNF^/zoc 

= F^Z/col. . . . . . . ( 16 ) 

The total intensity of the reflection is F'^Z/c, so its integral 
breadth in w is a. Thus the greater the chance of a mis- 
take the more the region of appreciable intensity of reflec- 
tion is spread out along the c* axis in reciprocal space. 

5. Having found the maximum intensity of reflection as 
a function of position in reciprocal space, that is G^{w) ^I{w) 
for zo = Oy let us try to evaluate I(zo) for an arbitrary value 
of w. In this case each layer will be reflecting ztzw out of 
phase with the preceding one, and the intensity of reflection 
as a function of w will be given by 

I(w) = G%zv) =A^+ 

where 

A = F^ + F2 cos znzv -h Fq cos + F^^ cos Sjtzo + • • *, 

B = F2 sin ZTtw -h Fq sin ^nzv + F4 sin Gjtw + ••*. 

By squaring and arranging terms as in equation (13) above 
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this becomes 

I{w) = Ff + + ••• 

+ 2[FiF2 cos 2,TCW + F2F3(cOS 2,7tW COS 

-f- sin 2.7ZW sin + • * *] 

+ cos /\Jtw + ^ 2 ^ 4(003 zuzw cos 67tw 

+ sin Z71W sin 6jcw) + • • •] 


= FI + FI + F| + ••• 

+ Z cos Z 71 W{FjF^ + F2F3 4- F3F4 + •••) 

+ z cos j\.nw(FjF2, + FJF'^ 4- FgFg 4 - •••) 

+ z cos 6 jtw\F^F^ 4- F2F5 4- FgFg + •••) 

4 * ***. ....... 

On replacing the products of the F’s by their mean values 
this becomes 


I{w) = NJq + z{N - i)/i cos Z 7 ZW + z{N - 2)/2 cos 4 ^^ 

4- z{N — 3)/3 cos 67 tw 4- . (18) 

where the /’s have the same meaning as before. This series 
also could be summed strictly for = F%i — aoc)”^, but 
for the present purpose it suffices to approximate the sum 
by an integral : 


roo 


where 


= zN 1 J.„i cos ZTtwm dm . 

J 0 

poo 

= aAT^F^i (i — zoif^ cos ZTVwm dm 

J 0 

poo 

= zNF^ I cos ZTtwm dm 

J 0 

y = — log ( I — 20c) *= 20c. 


(19) 

(20) 

(21) 


This integration may be done by parts : 


f e-^^' 

J 0 


—ax QQg — 


_g— ace ^Qg 


a 


=; sin hx 

a 


T-*fV 

Jo o 

'Jo 


sill bx dx 


cos bx dx. 



54 


X-RAY OPTICS 


^2 ! A2 poo ^ 

I cos hx dx — — , 

Jo a 


Q-ax _ a/(a^ + 


In the present case a = y, b = 'ztcw, so 

. zNF^y 

+ yzTtto)^ 


c(oC^ -t- TC^W^^} 

For zu o this reduces to equation (i6), as it should, 
course of liw) as a function of w is shown in figure 3. 
is very similar to that of sin^ TcNu/{nuYy 
shown in figure 2 of chapter III, except • |> # i 

that it does not oscillate. Both go to ' 



Fio. 3.— Variation of intensity as a 
function of zo for a layer structure 
■with, mistakes. 



Fig. 4. — Po 
of the reci 
cal lattice 
a layer S' 
ture with 
takes. 


zero as the inverse square of the variable for large va 
of the variable. The broadening of the reflections dil 
from that due to a small thickness of the crystal in 

those reflections with k odd 

^ected. A section of the reciprocal lattice is showr 
figure 4. 
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It is perhaps worth while to check that the total intensity 
of the reflection is actually ZF^/c, as this gives some assur- 
ance that the approximations used have not introduced 
serious errors. The total intensity is 

/ = I{w) dw = (xZF^/c)^^(ol^ + Jr2w2)-1 dw 

= (a.ZF^/c){Tca.)~^ arc tan (jtro/a)^ . 

At the upper limit tzw/ol has a large positive value, so that 
arc tan (jr^y/a) is practically tz/z, and at the lower limit it 
has a large negative value, so that arc tan (tczo/ol) is prac- 
tically ~ 7 t/z. The total value of the integral is therefore 

6. In the foregoing paragraphs the calculation of the 
intensity of reflec- 
tion from a layer 
structure containing 
mistakes has been 
carried through in 
detail for a particu- 
larly simple variety 
of mistake. The 
method adopted is, 
however, of general 
application in the 
calculation of the 
intensity of reflec- 
tion from a layer 
structure containing mistakes. If •**, ***, 

are the A and B parts of the structure factors of the unit 
cells in the various planes (equations (i)-(4) above), and c* 
is taken perpendicular to the layers, the contribution of 
the 72 ^^^ layer to the A part of the total amplitude of 
reflection is (figure 5) 

A^ cos Z7t{n ^ i)w — By^ sin 27 r(n — 



Fig. 5 . — Construction for the A and B 
parts of the amplitude reflected by 
the 71^^ layer. In this figure zrrnw 
should be Z'n(7i — 
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and the B part is 

sin 2.7t{n - i)w cos 27 i{n - i)w, 

and the intensity of reflection as a function of w is, the 
terms being arranged as before, 

I{p) == G\w) = A\w) + B\w) 

= (^A^ 4- A^cos'znzo + AqCosj{jzzo + ^40036712^; + 

- B^ sin 27CW - S3 sin ^tzzv - S4 sin 67tw - •••)^ 

+ (A2 sin 'ZTWO + A^ sin 471^ + Ai^ sin ^Ttw + 

+ Sn + SaCOSSTTZf; + S3COS47rzf? + S4 cos bTcet; + 

= (Af + S?) 4 - (Al 4 - S|) 4 - (Al 4- S|) 4- — 

+ 2 (Aj_A 2 4 - B1B2) COS ZTtW 

+ ziA^A^ 4 - S2S3)(cos ZTtW cos ^Tcw 4 - sin zTcw sin 477;^) 

+ 2(^3^ 4 4 - S3S4)(cos ^Ttw COS 6 nw + sin /[tzw sin (> 7 tw) 

4- *** 

+ z(AiAs + B^Bq) cos 4.Jtw - . x 

+ 2(^2-^4 + S2S4)(cos-27r2tJ COS 6jtw + Sin Z7ZW sin bTTW) 

+ 2(AaAs + ^gBsKcos 4yzw cos Snw + sin sin Sttzo) 

4 - **• 

+ ••• 

- 2 (^iS2 - ^a^i) sin . 

- zIAzBs - ^s^aXsin 47 cw cos znw - cos ^tiw sin ZTtw) 

- ziA^B^ - ^4jB3)(sin bitw cos ^nw - cos bjiw sin 477 :zy) 

• • • 

-* z(Aj^Bs - A3B1) sin 47 r 2 t; 

- 2 (^ 2^4 - ^4^2)(sin 6Ttzo cos ZTtw - COS 6tiw sin znw) 

- z{aIb^ - ^gjBaXsin Snw cos ^w - cos Sttz^; sin 477:7^) 

■ 4 • • 

• 4 • 

= {A\ 4- 4- {A 14- ^1) + {Al 4- ^1) 4- 

4 - 2 cos 27 rZ 7 j[(/ 4 i ^2 + B^B,^ -I- (A2A3 4 - B2B2) 

4 - (A^A^ 4 - B^B^) + — ] 

4- 2 cos ^Ttw\J{Aj^A3 4 - B^B^ 4- (A2A 4- B2B 4) 

4- (^3^5 + -^3^5) + 

4- 2 cos 6ftw[{A^A^ + B^B^) + (AqA^ 4 - 
+ (A^Aa 4 - ^3-60) 4 - ***] 
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— 2 sin 2 , 7 czu[{A^B 2 — + {A^B^ — -^3^2) 

+ (^ 3-^4 - A^B^ + -••] 

— 2 sin 47 i:zy[(^i-B 3 — A^B^ + {A^B^^ — A^B^ 

+ (^ 3^5 - A^B^) + -] 

•**, C^s) 

This equation contains terms of two kinds, those of the 
form AA' + BB' , and those of the form AB' — A'B, where 
A and B belong to one layer and A' and B' to another 
layer. If is the mean value of terms of the first kind 
when the layers are separated by m interlayer spaces, and 
the mean value of terms of the second kind, equation 
(25) becomes 

I{w) = N J q+^(N cos Z 7 tW+z{N — ‘Z)j2 cos /[TIW'' 

-H ••• +2 Jjsf-1 cos ZJc{N — l)w 1 / ^\ 

-~2(iV— i)i^i sin 2 jr 52 ?— 2(iV— 2)i^2 4^^ j ^ 

— ••• — zK]^_-i SAX). Z 7 t{N— 1)^0. J 

This is the generalized form of equation (18), and is valid 
in every case when the number of layers is great enough 
to allow averaging the values of AA' + BB' and AB' — A'B. 
If in addition the chance of a mistake is large compared 
with I /AT, but small compared with j/z, the sums may be 
approximated by integrals, so that 

rN 

I(zo) =21 (N — cos zTtwm dm 

Jo ^ 

— 2 I (iV* — sin ZTtzom dm. . (27) 

J 0 

The correct upper limit is iV = Z/c, but with the chance 
of a mistake within the range stated it can just as well be 
made infinity, and {N — m) can be replaced by iV. In 
many cases therefore equation (26) may be replaced by 

pQO 

= zN I COS ZTZwm — sin ZTtwrn) dm^ (28) 

J 0 

but it must be remembered that this is an approximation, 
and if the chance of a mistake is very large or very small 
equation (26) must be used as it stands, or approximated 
5 
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in a way appropriate to the particular problem. From the 
definitions of and it is obvious that J _ ^ == J and 
K = — K so that equation (28) may be written more 

symmetrically as 

r+ °° 

I{w) = N I (Jm COS 27 cwm - sin znwm) dm 


or 

I{zo) = cos 'ZTVwni — sin 2 . 7 vwm) dm . (29) 

for short. In very many cases, as in the example used in 
paragraphs 4 and 5» i® zero and the terms involving it 
vanish. In most of the problems to be considered in the 
following chapter J^n has the form ^.s in equation 

(21) above, and equation (28) becomes 


I(w) - 


+ {2.7tzo)^ 


• (30) 


The maximum value of this is zNJq/y^ and the total 
intensity is NJq, so that the integral breadth in zi; is y/2. 

7. The general effect of mistakes in layer structures will 
now be obvious. Reflections whose structure amplitude 
is not affected by the mistake in question (for example, all 
reflections with k even in the example considered above) 
will be unaffected, and the spots or lines on the X-ray 
photograph will be as sharp as if the crystal were perfect. 
On the other hand, reflections whose structure amplitude 
is affected by the mistake will be drawn out in reciprocal 
space in a direction perpendicular to the layers. The 
amount of the spread will depend on the frequency of the 
mistakes and the amount of the change in the structure 
amplitude ; in the simple example considered in para- 
graphs 4 and 5 ^11 reflections with k odd were drawn out 
into streaks in reciprocal space of length of the order of 
c^/oc. The corresponding spots on the X-ray photograph 
will in general be drawn out also, or the lines in the powder 
photograph broadened, but the magnitude of the effect 
depends on the angle that the streak makes with the 
surface of the sphere of reflection, and is most readily 
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discussed in connection with particular examples, as is done 
in the following chapter. If is not zero the reflection 
is displaced as well as broadened. This is readily seen 
from any of equations (26) to (29), since /(- w) is not equal 
to I{wy unless the terms involving sines vanish. One or 
two examples of displacement are discussed in chapter VII. 

8. Besides layer structures with mistakes, which may be 
described as mistakes in one dimension, examples are known 
of substances exhibiting mistakes in three dimensions. One 
of these is an alloy of gold and copper, in the proportion of 
one gold atom to three 
copper atoms. The nature 
of the mistakes may be 
understood from figure 6, 
which represents one unit 
cell of the crystal. The 
cell is cubic and contains 
four atoms, distinguished 
by the letters J 5 , C, D. 

(At first sight there appear 
to be twelve atoms, but 
only one-eighth of each' 
atom marked D is really inside the cell, the remaining 
seven-eighths being in neighbouring cells. Similarly, 
only one-half of each atom marked A, or C is inside 
the cell.) On the average, one of these four atoms will 
be gold and three will be copper. At high temperatures 
the gold and copper atoms are mixed at random, and 
powder photographs show only the reflections corre- 
sponding to a face-centred cubic structure. At lower 
temperatures the gold atoms, behaving as if they were 
larger than the copper atoms, avoid one another, and in 
an ideal crystal would all settle into position Aj or all into 
position jB, or all into position C, or all into position Z>, 
the three positions not occupied by gold atoms being 
occupied by copper atoms. Thus at high temperatures 
any position, say A, is occupied on the average by one- 
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quarter of a gold, atom and three-quarters of a copper 
ato'm, and at low temperatures either always by gold or 
always by copper. There are therefore at low temperatures 
four distinct types of unit cell : 

jiy with gold atoms in the centre of the face perpendicular 
to the a axis, 

JB, with gold atoms in the centre of the face perpendicular 
to the b axis, 

C, with gold atoms in the centre of the face perpendicular 
to the c axis, and 

Z>, with gold atoms at the cube corners. 

When the gold atoms have sorted themselves out in one 
of these ways throughout a crystal it is said to be ‘ ordered \ 
and a large number of extra X-ray reflections appear, cor- 
responding to the decrease in symmetry from statistically 
face-centred cubic to simple cubic. If the ordering is 
complete, each gold atom is entirely separated from con- 
tact with others by a screen of copper atoms. A com- 
paratively small increase in temperature is able to overcome 
the gold atoms’ dislike of one another’s company, so that it 
would not be surprising if within an ordered single crystal 
mistakes occurred such that the type of order varied from 
place to place. This is found to be the case, and by proper 
heat treatment it is easy to put the alloy into a state in 
which the reflections common to the ordered and disordered 
state are sharp, but those characteristics of the ordered state 
only are more or less diffuse. One or two other substances 
exhibiting mistakes of this type are known, but they are 
much rarer than layer structures with mistakes, perhaps 
because it is only possible to build up three-dimensional 
mistakes with exactly four types of unit cell (Bragg, 1940). 

9 * The intensity of reflection from such a structure as 
a function of position in reciprocal space could be treated 
by methods similar to those of paragraph 6 above, but the 
double series there, already sufficiently clumsy, are replaced 
by sixfold senes, and become almost unmanageable. Con- 
sideration of this problem will therefore be postponed to 
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chapter VII, where more advanced mathematical methods 
are introduced. Line broadening in powder photographs 
may, however, be treated here by essentially the same 
argument as was applied previously in connection with 
particle size broadening (chapter IV). To calculate the 
integral breadth of the hkl reflection, we calculate the in- 
tensity scattered with *9 between *5 and + dS by a single 
line of unit cells parallel to the radius vector joining the 
point h, k, I of the reciprocal lattice to the origin, and sum 
for all parallel lines. If paragraph 6 is re-read, it will be 
seen that the calculation of the intensity scattered by a 
single line of cells goes through in exactly the same way 
as the calculation of I{w), except that and are the 
mean values of A A' + BB' and AB' — A' B for cells separ- 
ated by m repeat distances parallel to the line, instead of 
perpendicular to the layers. With axes chosen as in 
paragraph 3, chapter IV, one may rewrite equation (27) as 

rT 

H(u) du {2,/ a)\ {T — cos {zntu/ a) dt du 

J 0 


rT 

— (2/<2)l (T — t)K^(t) sin (zTctu/a) dt du . (31) 

J 0 

where a is the repeat distance in the direction of the line, 
T = Na is the length of the line, t = mUy and u is written 
in place of w to conform to the notation of the previous 
chapter. The interpretation of J{f) and Kif) is fairly 
obvious ; they are the mean values of A A' + BB' and 
AB' — A' B for cells separated by a distance t along the 
line. Equation (31) is analogous to equation (3) of 
chapter IV, and ^(6') is obtained from it in the same way : 


I^S) = (a/U)^^H{u) dy dz 

■ [T — cos {■ZTztu/d) dt dy dz 

— t)IC(t) sin {Z7itii/d) dt dy dz 


(32) 
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By equation (15), chapter IV, the volume common to a 
crystal and its ‘ ghost ’ shifted a distance t in the hkl 
direction is 

V{t) = JJ(T - t) dy dz, 
so 

= (2/t/) [^V{t)J{t) cos (2,7Ztu/a) dt 

J 0 



sin {z 7 ztu/d) dt 



V {£)!{€) cos ZTcst dt 


— { 2 /U){ V{t)K{t) s\ll2,7ist dt . (33) 

J 0 

where the integration extends to a large enough value of t 
for V(t) to vanish, and ^ — uj a is the actual distance by 
which S misses the point A, k, I of the reciprocal lattice. 
It will be remembered that u is the fraction of the distance 
I /a between the reciprocal lattice points by which S misses 
the point of the reciprocal lattice. Since J{t), V(t) and 
cos 27tst are unaffected by a change in the sign of ^, and 
J^(t) and sin zjtst both change sign, equation (33) may be 
written 


I(S) = U-^[fr(t)J(t) cos 27 tSt dt 

— iV{t)K{t) sin 27 zst dt] . (34) 

where the limits of integration are the negative and positive 
values of t for which V{f) vanishes. This equation is very 
similar to equation (21) of chapter IV, and in fact reduces 
to it if there are no mistakes, for then J{t) is constant and 
K{t) is zero. 

Ordinarily J{t) will have the form 

Jif) = Joexp(-y ji| ), 

where y is a function of h, k and /, and Kif) will be zero. 
If mistakes are sufficiently frequent to cause line broaden- 
/(^) will become inappreciable before V{t) departs 
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greatly from its value for i = o. Equation (34) then 
becomes 

p+co 

7 (*S) = JqVU-^ I exp {—y\t\ ) cos znst dt 

J —00 

_ JoV zy 
U + {zTts)^ 

The maximum value of this is 1 ^ zJ^V/Uy, so that the 
apparent particle size is 

e=IJI = z/y,. . . . (36) 

and the integral breadth is 

P - X/e cos 6 = Xy/z cos d. . . (37) 

When K{t) is not zero one may expect displacement as 
well as broadening of the lines. 



REFERENCES 

W. L. Bragg, 1933. The Crystalline State. London : Bell and 
Sons. 

1940. Proc. Phys. Soc.y 52, 105. 

S. B. Hendricks and E. Teller, 1942. y. Chem. Phys.y 10, 147. 
Internationale Tabellen sur Bestimrnung von Kristallstruktureny 
1935. Berlin : Borntrager. 

I.r. Landau, 1937. Phys. Z, d. Sow. Uti.y 12, 579. 

1 . M. Lifschitz, 1937. Phys. Zeit. d. Sow. Un.y 12, 623. 

1939- J' Pxp- Theor. Phys. USSR.y 500. 

K. Lonsdale, 1936- Simplified Structure Factor a?id Electron 
Density Formulae for the 230 Space Groups of Mathematical 
Crystallography . London : Bell and Sons. 

A. J. C. Wilson, 1942. Proc. Roy. Soc.y A, 180, 277. 

1943- Proc. Roy. Soc.y A, 181, 360. 

(References to experimental work on substances mentioned 
above will be found in the next chapter.) 



CHAPTER VI 


MISTAKES IN LAYER STRUCTURES : SILICATES, 

COBALT, GRAPHITE 

MISTAKES IN THREE DIMENSIONS : POWDER 

PHOTOGRAPHS OF AuCug 

I. The mistake specially considered in paragraph 3 o£ 
the preceding chapter does not seem to occur, but a 
similar one with a shift of one-third or two-thirds of 
the b axis is found in some layer silicates (Hendricks, 
1940). Suppose that the three types of layer are A, 
and C, and that the crystal starts off AAA*‘% but that 
at each added layer th^re is a chance a of a mistake 
occurring, and the sequence continuing BBB ‘* * or CCC* • *. 
Suppose, further, that either mistake is equally likely. 
Let the probability of the added layer being A be 
Any of being B be Bn, and of being C be 

Then An+i = (i - cx)An + i-cxBn + JaC^, . (i) 

the factors of one-half coming in because, if the added 
layer is R or C and a mistake occurs, the {n + will 

be A only half the time. Since An + Bn -H = i, 
equation i may be written 

■^w+l “ (t 

This may be solved in the same way as equation 6 of 
the preceding chapter, the solution being 

i + 0(i - fa)" . . (3) 

where Q is an arbitrary constant. Similar equations can 
be set up for Bn and On, with the same solution. Choosing 
the constants so that Aq i, Bq = Cq = o, gives 

A. = i[i + 2(1 - foe)-],'] 

A == i[i - (i - foef ], V . . (4) 

c. ^ i[i - (I - fa)-]. J 

If the structure amplitude of the A plane is F, that for the 

64 
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lave an ^ part of F cos i^vzk/'^ and a JB part 
/3), and that for the C plane will have an 
cos (4:7r^/3) and a JB part of F sin (4^-fe/3). 
mean value of the product JLJL' + FJB' is 
the fact that ^ = o, so that is the mean 
and is given by 

cos (ajr^/s) + cos (4^^/3)] 

+■ cos (jZTtk/^) + COS (47rA/3)} 

-I- (i — ^oc )”^{3 — cos {p.Ttk/'^') — cos (4:7r^/3)}] 

^ = 3 X an integer, 

^ = (3 X an integer) =b i- 

4uation (5) it was tacitly assumed that one of 
ncerned was an plane, but the same result 
' a B plane or a C plane is assumed. It is 
hat is zero : 

+ B^ sin ( 2 , 7 tk/ 2 ) + sin ( 47 r^/ 3 )] 



and sin (43 t^/ 3) = — sin ( 27 t^/ 3 ). 

= F^ the reflections with k/'^ integral are 
ut those with ^/3 fractional are drawn out 
n. reciprocal space parallel to the c'^ axis, the 
intensity along the streaks being given by 
, chapter V : 


lizo) 


'zlSfF’^y 

laA^-F^oc 

9CC^ + 



/s. On a rotating-crystal photograph these 
streaks along curves in the photograph corre- 
Lines parallel to the axis, but on a powder 
they will be somew^hat fore-shortened, and 
:h C>Q is perpendicular to will be practically 
For example, in an orthogonal crystal the 
= o will be sharp, as well as those with k 
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divisible by three. A more extended discussion of the fore- 
shortening is given below, in connection with line-broaden- 
ing in photographs of cobalt, and of the lines with Z = o, 
in connection with graphite. Hendricks (1940) has made 
an extensive investigation of silicates showing mistakes of 
this type, but has not attempted a quantitative comparison 
with theory. 

The other mistake possible with the silicate minerals is 
an. increase in the interlayer spacing at random intervals. 
If the mistake consists of an increase of the spacing by 
a fraction g, each mistake will cause an increase of the phase 
of all succeeding layers of 'ZTtlg. The mean value of the 
product AA' + BB' for layers separated by m interlayer 
spaces will therefore be 

~ multiplied by the chance of no mistake in m layers 
+ cos iZTtlg times the chance of one mistake in 
m layers 

+ F^ cos ^Tzlg times the chance of two mistakes in 
m layers 

• • • 

+ F^ cos 'ZTcmlg times the chance of m mistakes in 
m layers 


= + ma(i — cos ZTilg 

4- - - cos ^Tzlg + ••• + cos 

3! 



In the problems so far considered the mean value of 
AB' — A'B, has been zero, a considerable simplification 
being thus produced. In the present case, however, 
does not vanish, but is given by an expression like that for 
except that the cosines are replaced by sines, so that 

Kjii = sin 2.Tclg 


fjzCffz — 

+ i(x^(i — sin 4zc/p + •** + a’” sin ^nnilg 

3! __ 
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(The term in (i ~ is missing in the expression for 
because it is multiplied by the sine of zero, and so 
vanishes, whereas in it is multiplied by the cosine of 
zero, which is unity.) With expressions as complicated as 
these for and it is not possible to evaluate equation 
(z8), chapter V, simply by elementary means, and recourse 
must be had to methods involving complex quantities. It 
may be said, however, that different from zero means 
that zo = o is no longer the position of maximum intensity, 
so that we may expect a change in the positions of the 
reflections as well as a broadening. 

2. In hexagonal cobalt the mistake is a shift of one-third 
in the a direction and two-thirds in the b direction, or else 
twice this (paragraph 2 , chapter V). If the structure 
amplitude of an A plane is F, that for a B plane has an 


2k\ ^ fh—JC 


A part of F cos znl — i- — ) = F cos ztt 

3/ V 3 


I, and a B 


part of F sin Z7t{ ~ = jPsin ~ ]. while that 


O' 


of a C plane has an A part of F cos ztz 


'h - K 


and a B part 


of — F sin Z7z{ 


fh - K 


\ 3 


I. The problem differs from that 


just treated in that the layer and the {n - 1 - must be 
different ; the mistake consists of the (w + 2 )^^^ being unlike 
the Thus if the first plane is A and the second B^ the 

third ought to be A, but there is a chance a that it is C. 
Let An be the chance that the added plane is like the 
first. The probability that the (zz + 2 )^'^ is like the (w + 
is zero, that it is like the is (i — oc), and that it is like 
neither is a, so 


O X + Oc(l An "^ri-l-l)) 

or An^.^ + - (l - ^CL)An = OC. . (7) 

This is a nastier equation than either (6), chapter V, or 
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(2), above, but it can be solved in the same way. Let 

= -P + 

P + — (1 — 2a)P - (r - 2oC)Qp”' = «, 


P = i/ 3 > 


-. 7^+2 


+ ap 


«■+! 


(i — 2 a)p^ = o, 


— oc =t: '\/a^ H- 4 — 8 a 

P . 

2 



It will be found that the positive sign leads to reflections 
that, with the ordinary hexagonal indices, have I even, 
and the negative sign to reflections that have I odd. Let 


Then 


— a + A/4 — 8a + a^ 

Pe — — 5 

2 

— a — A/4 — 8a + a^ 

Po *“ ~ — • 

2 

= h + QeP^ + QoPJ^ 


(9) 

(10) 


where Qg and are arbitrary constants. They must be 
chosen so that Aq = x, = o ; the values required are 
readily found to be 


^Po')/ 3 (.Po Pe),\ 

go = (X + ^Pe)/3{pe - Po)‘j 
For c/.Z' « a <<( I these reduce to 

Pe = I - fa N exp (- 3 a/ 2 ), Qe = i/6,l 
Po = - X + *=. - exp (- a/2), ge = X/2.J 

The values of and C^, the probabilities that the 
added plane will be B or O, are given by equation (10) 
with the constants chosen to make Bq = Cq = o, 
-Pi = C*! = 1/2, since the first added plane is equally likely 
to be B or C. It is found that 


-Pn = i + -P^P/' + ^oPo^ 
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where 

“ (l + 2pJ/6(p^- Po) - l/l2,\ 

“ (l + 2pg)/6(^^ Pg) == '^/A' J 

The mean value of AA' + JBB' is then 



(13) 


T = F^A + F^B 



4- F^C^^ cos Z7t 



= F^lUi H- zC) + (Q, 4- zCR,)pr 

+ {Qo + • * • (14) 


where C has been written for cos Z7t 



For h — k 


divisible by three this reduces to F^ ; for h — k not divisible 
by three it reduces to 


Jm. = 

or exp (- 3am/2) 4- | cos mTt exp (- ocw/2)] (15) 

approximately. The value of is readily found to be 
zero. Reflections with h k divisible by three are normal, 
since for them — F^ ; the line-profiles for the other 
reflections are given by 

pOO 

I{tJo) = zA/’i cos ZTtwm dm 

J 0 

zNF^ I {J exp ( — 3am/z) 

J 0 

4- ^ exp (— (x,m/z') cos tHTt} cos Z 7 zw 7 n dm 


= zNFV - 


4/30^' 

I VT 


30^ 


4- {ZTtw)^ 


a 

z 


8/oc 



4- [zn{zo 4- J)]^ 


a 

z 


8/oc 



>- 


2 4 - [ZTC{W - -i)]2 


(16) 
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C^tie second integration is done by expressing the product 
of two cosines as the sum of two cosines.) This is a rather 
more complex expression than we have obtained previously, 
but th.e interpretation is straightforward. The points of 
the reciprocal lattice have broad weak peaks of integral 
breadth a = 3 a/4 and intensity while halfway 

between them, there are sharper stronger peaks of integral 
breadth cr == a/4 and intensity 3iVP'^/4. The reason for 
this apparent halving of the reciprocal lattice is that we 
have been taking a unit cell as consisting of a single plane, 
whereas the normal unit cell of hexagonal cobalt consists 
of two planes. The reciprocal lattice, wdth the ordinary 

values of /, is illustrated in 
f figure I. The points with 

^ h — k divisible by three and 
I I even are sharp, and with I 
f odd are missing. The 

I points with h — k not divis- 
^ ible by three and I even are 
I drawn out into long streaks 
f parallel to the c* axis ; those 
with I odd are also drawn 
out, but only one-third as 
much. The integral breadths 
referred to the ordinary values of I are twice those given 
above, so that g = sa/z for reflections with I even, and 
<y = oc/3 for reflections with I odd. 

In the case of cobalt the strict calculation of the intensity 
has been carried out, equation (26), chapter V, for the line- 
profile and equations (9), (ii) and (13) for the p’s and Q’s 
being used. For a small compared with c/Z the reciprocal 
lattice for normal hexagonal close-packing is obtained, for 
a large compared with c/Z but small compared with one- 
half it is as just described, and for a t/a the fuzzy 
regions practically disappear. For a a little greater than 
one-half the p’s become imaginary, but I{w) can still be 
interpreted, fuzzy regions appearing in places correspond- 


I 




I 




i! 


■4 




% 


I 




I 


Fig, I . Portion of the reci] 

rocal lattice for cobalt. 
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ing to the points of the reciprocal lattice for cubic close- 
packing. As oc approaches unity these regions coalesce 
into sharp points, the final state being the superposition of 
the reciprocal lattices for the two orientations of the cubic 
crystal corresponding to the sequences of planes 
ABCABC*'' and ACBACB^-^. This is as it should be, 
since an hexagonal crystal with a mistake at every plane 
(oc = i) would be cubic close-packed, but the beginnings 
AB and AC are equally likely. 

A . careful experimental study of the broadening of the 
X-ray reflections from imperfect hexagonal cobalt has been 
made by Edwards and Lipson (1942), both single-crystal 
and powder photographs being used. The single-crystal 
photographs showed in a striking manner the extension 
of the regions of reflection along lines parallel to the c* 
axis, and the fact that reflections with I odd were much 
less affected than those with I even. They were, however, 
unsuitable for quantitative measurements, and for this pur- 
pose a powder photograph taken with crystal- reflected 
nickel radiation was used. (Ci*ystal-reflected radiation 
is practically monochromatic, 
consisting almost entirely of 
one doublet, and so gives 
photographs with less back- 
ground than merely filtered 
radiation does.) The integral 
breadths were derived by the 
method of Jones (1938)- The 
integral breadths are given in 
the second column of table I. 

They depend, however, on 
the actual value of Z, as well 
as on a and the parity (odd- 
ness or evenness) of Z. This is easily seen from figure 2. 
A ‘ spike ’ of length in reciprocal space <5 corresponds 
for Z = Zi to a range of sin 6 proportional to B^Cj^y whereas 
for Z = Z2 it corresponds to the greater range proportional 



Fig. s. — Dependence of the 
breadth of a reflection on Z. 
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to B^C2‘ The relation between and I is most easily got 
from the equation 

sin.^0 = + hk + • (i?) 

which relates the angle of reflection to the indices of 
reflection for an hexagonal crystal. Differentiation gives 

2Sin0cos0— = — - • (18) 


If G is the integral breadth in /, or which amounts to the 
same thing, 


P — 2(7 


dQ 


G 


dV 

== = 


sin 20 
“1 


/ 5 . 


• (19) 


TABLE I 

Calculation of a, the probability of a mistake per added plane of 

atomSj in a specimen of Cohalt 


Line 


cr 

a 

a corr. 

loi I 

0*0092 

0*043 

0*086 

o*o 88 

1012 

0*0320 

0*089 

0*059 

0*073 

1013 

0*0256 

0*051 

0*102 

0*1 1 1 

2021 

0*0076 

0*044 

0*088 

0*093 

2022 

0*0406 

0*105 

0*070 

o*o 88 

1014 

0*0844 

0*099 

0*066 

0*1 19 

2023 

0*0592 

0*042 

0*084 

0093 


It follows that p = o for I = o, and that for a given value 
of 6 it is smaller the smaller the value of /. The experi- 
mental values of a, derived from those of /5 by means of 
equation 19, are given in column three of the table. It 
will be seen that they fall into two groups, those with 
G> — '0*045 odd), and those with — 'O*io (/ even). 

The theory given above predicts g = 3a/2 for the lines 
with I even, and g = ol/z for the lines with / odd. The 
value of cx, the chance of a mistake per added plane, is 
thus about 0*090 if calculated from the lines with I odd, 
and about 0*067 if calculated from the lines with I even. 
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This is not entirely satisfactory agreement. Though the 
value of a is not very small compared with one-half, 
errors due to the various approximations made (for 
example, the use of equations (12) instead of (9) and (10) ) 
do not seem sufficient to account for the discrepancy 
directly. A more probable explanation is limitation of 
the range of photometry. The theoretical profile of I(zjo) 
is shown in figure 3. In calculating the integral breadth 
as cr = // dw/I(p) it is tacitly assumed that it is possible 
to photometer a sufficient length on the film on each 
side of the line to get practically all the intensity. In 
the present case the region that it was possible to photo- 



Fig. 3 . — Correction of the intensity for limitation of the range 

of photometry. 

meter corresponded approximately to the range from 
= — 0*15 to et?=+o*i5, so that the intensity corre- 
sponding to the shaded portion of the figure was lost. 
A correction can be calculated (Wilson, 1942), and it is 
found that lines with I even are aflFected more than the 
lines with I odd. The corrected values of a are given 
in the fifth column of the table, and it will be seen that 
the discrepancy has practically disappeared. The mean 
of the corrected values of a is 0*096 ± o-ooS from the odd 
lines and 0*093 =fc 0*014 from the even lines. It seems, 
therefore, that the assumption of mistakes in the structure, 
distributed at random along the hexagonal axis, explains 
satisfactorily the broadening of the lines and accounts 
6 
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quantitatively for the difference between the lines with 
different parity. 

The correction to the breadths may be calculated as 
follows. If the integral breadth in w is c. 


I{w) = 





(equation 30, chapter V, with y replaced by 2cr). If the 
range of photometry is from w ^ — x to w + x, the 
observed maximum intensity will be 

Iq * /(o) — I{x) 

^ + Ti'^x^) . . ( 21 ) 

instead of = /(o) = NJ^/g. . . . (22) 

The observed total intensity will be 


r = 



dw — zxl{x) 




uz 




arc tan {tzx / a) 


Tcxa 


G^ -h (jCX)'^^ 


p+00 

instead of 7=1 I{w) dw = NJq, 

J —00 

The observed integral breadth will therefore be 


2cr 


= — [ (i + g^^/tz^x^) arc tan {tzx/g) — G/nx"], 


(23) 

(24) 

(^ 5 ) 


instead of cr as it ought to be. For x =0*15, g' — 0*045 
(mean for I odd) a little trial and error gives g = 0*051, 
and for g' = 0*098 (mean for I even) it gives a = 0*140. 
The corresponding values of a are 0*102 and 0*093, in fair 
agreement with one another and with the values found by 
AVilson by correcting each line individually. 

3 * -Another substance exhibiting mistakes is graphite, 
and it has been the subject of numerous investigations in 
connection with lamp black, coke, arc carbons and other 
commercially important products. The ideal graphite 
structure consists of layers like those shown in figure 4, 
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alternate layers taking the positions A and B. Well- 
crystallized graphite usually consists chiefly of this struc- 
ture, but contains 
also an appreciable 
fraction of another 
structure in which 
the layer sequence 
is ABCABC — 

(Lipson and 
Stokes, 194Z). The 
relation between 
these structures is 
similar to that 
existing between 
hexagonal and 
cubic cobalt, but intermediate structures do not seem to 
have been reported, and the remainder of this discussion 
will be confined to disorder in the ideal form. The binding 
between adjacent layers is weak, so that they are readily 
displaced or rotated relatively to one another by any 
amount. This behaviour contrasts with that of cobalt and 
the layer silicates, which allow only a small number of dif- 
ferent displacements. In addition, the crystals are small, 
so that it is necessary to consider ‘ particle size ’ broadening 
as well as ‘ mistake ’ broadening. 

The structure of ordinary, non-ideal, graphite may be 
pictured as a sequence of planes ABAB--- for a certain 
distance, but at each added plane there is a chance of a 
random translation and rotation, destroying completely the 
regular phase relation between atoms separated by the dis- 
placement. This will not affect the 00/ reflections, so that 
for them the intensity will depend on w through the 
ordinary factor 



I{w) 




sin^ TtNw 




(26) 


X-RAY OPTICS 


76 

where N is the number of unit cells in the crystal in the 
c direction (chapter III), but all other reflections will be 
drawn out parallel to the c* axis by an amount depending 
on a. The calculation of is complicated slightly by the 
fact that a unit cell consists of two layers. If a is not too 
large, however, it will be near enough to take the chance 
of a mistake per unit cell as 2a, and treat the crystal as 
consisting of N layers of unit cells, instead of zN actual 
layers. Then the mean value of A A' + BB' is if no mis- 
take intervenes between the layers of unit cells concerned,* 
and zero if a mistake does intervene, since there is no 
regular phase relation between layers separated by mis- 
takes. The value of is therefore multiplied by the 
chance of there being no mistake in m added layers of 
unit cells, so 

= F\i - zoc)^ '== F^ exp (- zcam). . (27) 

The intensity of reflection for reflections other than the 
oo/’s will therefore depend on w through the factor 


, NolF^ 

zo) = . 

+ (tzzv)^ 

(equation (24), chapter V). For a not too large a specimen 
consisting of such crystals will give ool powder lines as 
sharp as the length of the crystal in the c direction permits, 
and hkl lines with breadths depending on a and I through 
a relation similar to equation (19) : 

XH 

sin z^^ 





the integral breadth in w of (28) being a. 

The reciprocal lattice just described, consisting of rela- 
tively sharp points on the axis, and all other reflections 
represented by streaks, suffices to treat powder photo- 
graphs, but it is not strictly correct. The regions separated 
by mistakes are not only translated with respect to one 
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another, but are also rotated by random amounts. This 
corresponds to a rotation of the reciprocal lattice about the 
c* axis, so that the reflections, other than the ooZ’s, will be 
turned into cylinders with centres on the c* axis, the inten- 
sity being uniform around the circumference of the 
cylinders, but depending on to according to equation (28). 
A single-crystal photograph of graphite containing mistakes 
of the sort considered here would be expected to show, 
therefore, broad bands in directions perpendicular to the 

axis, rather than narrow streaks parallel to it. 

For large values of a the approximations leading to 
equations (28) and (29) are no longer valid, and further 
consideration is needed. It does not seem worth while to 
treat intermediate cases, but the final state of a shift at 
every added plane (a = i) is of importance. The ooZ 
reflections will be, as before, relatively sharp points depend- 
ing on to according to equation (26), but the other reflec- 
tions will be joined up into continuous lines parallel to 
c, and of intensity uniform in to except for the falling off 
of the atomic scattering factor with increasing sin 6 /X. The 
crystal behaves essentially as if it were single for the ooZ 
reflections, giving points as sharp as its size allows, but as 
an assemblage of non-co-operating planes for all other 
reflections, giving lines in reciprocal space perpendicular 
to the planes (chapter II, paragraph 4). Crystals as faulty 
as this will not be large enough to give single- crystal photo- 
graphs, so that the rotation of the reciprocal lattice, pro- 
ducing cylinders instead of lines, need not be considered. 

A powder photograph of a specimen consisting of such 
crystals will show lines of two types. Lines of the first 
type correspond to the ooZ reflections, and are essentially 
similar to the broadened lines considered in chapter IV. 
The apparent particle size may be got from equation (26), 
but it is perhaps instructive to use the method of chapter 
IV. For this purpose the crystal may be considered 
a cylinder (not necessarily circular) of length Z, so that 
F'(Z), for a shift parallel to the c axis, is P{Z — |^| ), where 
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P is the cross-section of the cylinder, so that the apparent 
particle size is 

e = V-^^V(t) dt 


= P{Z -\t\)dt 


z. ...... ( 3 ^) 

Thus the apparent particle size calculated from the integral 
breadth gives the true length in the c direction. For 
comparison with reflections of the second type, to which 
the concept of integral breadth is not applicable, it is inter- 
esting to work out the half-width of the 00/ reflections. 
For these the line profile is given by 


I{S) 


sin 2 TtZs 
{ns)^ 


(31) 


where proportionality factors have been omitted, 
S = z sin 6 /X is the distance from the origin in reciprocal 
space, and s is the amount by which S differs from its value 
*So at the point ooZ. The intensity has the value Z"^ for 
5=0, and falls to half this for the value of s satisfying 


or 


sin^ TcZs 
{utsY 



5 = zh 0-442/Z' 


(32) 


approximately. Thus the half breadth in s is o- 88 /.Z', and 
the half breadth in zO, the quantity corresponding to 
integral breadth, is A/cos 6 times that in 5, or 

0-8851 




Z cos 0 


(33) 


The apparent particle size calculated from the half 
breadth is 

s' = X/fi' co&e = I-X-^Z, . . (34) 

and is thus 13 per cent large. 

Lines of the second type consist of an amalgamation of 
all hkl reflections having the same values of h and k, and 
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differ from any so far encountered. The calculation of the 
line profiles involves mathematical functions beyond the 
scope of this book, so we must content ourselves with 
a qualitative treatment. As we have seen, through every 
lattice point (A, k) in the plane I — o there passes a rod- 
like region of reflection. The intensity on the film with 
6 between 0 and 6 + dd will depend on the amount of this 
rod that lies between concentric spheres of radii S and 
S + dS, where dS = 2 , cos d dd/X (figure 5). The radius 
of the rod is of the order of i/2/>, where p is the square 
root of the area of a layer. Let s he the amount by which 



Fig. 5. — Calculation of the line profile for hkl reflections. 


S differs from its value Sq at the point hko. The intensity 
due to the rod will be zero until S is approximately 
Sq — 1/2^, when it will increase rapidly until S = Sq. 
For S^Sq the volume of the rod intercepted by the spheres 
continues to increase until S is approximately Sq + 
Beyond this point it decreases again, but less rapidly than 
it initially increased. For S > Sq the intensity on the film 
will at first continue to increase with the increasing volume 
intercepted by the spheres, but will not do so right up to 
Sq + 1/2^, because the intensity in the rod is greatest in 
the middle and falls off at the outside. The maximum on 
the film will therefore lie somewhere between the value of 
d corresponding to S ^ Sq and that corresponding to 
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or at about *So + i/ 4 ^> say. The peak of 
the reflection is therefore displaced, from the position where 
it would be expected by 

Ad = (A/a cos 0 )(i/ 4 p) = 0-125 X/p cosS. . (35) 

This displacement was discovered experimentally before 
the peculiar nature of the reflections was realized, and was 
interpreted as a contraction of the lattice in small crystals 
of graphite. 

The approximate line profile as a function of ^ is shown 
in figure 6 . The concept of integral breadth is not useful 
for dealing with such reflections, as the intensity never 
falls quite to zero with increasing and it is necessary to 



use the half breadth as a measure of the broadness of the 
reflection. The intensity reaches its maximum for 
^ '^/ 4 -p, but from the volumes intercepted by the 

spheres in the two positions it is pretty clear that the inten- 
sity for = o is greater than half that for ^ = i / 4 p. At 
a guess we might say that the intensity is half the maximum 
at the same distance to the other side of 5 = 0, or at 
s = — i/^p. The distance on the low-angle side between 
the maximum and the half value is thus i/2j>. On the 
high-angle side there is nothing to aid in the estimation 
of the distance it is necessary to go to reach the half value, 
but as the rise in intensity is steeper than the fall it will 
certainly be more than on the low-angle side, and as a guess 
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we might double it and say i /p. The total half breadth in 
s of the reflection would then be 3/2^, and the half breadth 
in zd would be 

/ 3 ' = (A/cos 6)(3/2^>) = i‘5Vi> COS0, . (36) 

and the apparent particle size calculated from the half 
breadth would be 

e' = X/P' cos 0 = o-byp. . . (37) 

If our estimates are anywhere near right, therefore, the 
numerical factor is considerably different from the 1-13 
derived above for the 00/ reflections. Warren (1942) has 
calculated more exact values, assuming that the intensity 
falls off exponentially with the square of the distance from 
the centre of the rod, and obtaining 

E — o*57P • • • • ( 3 ^) 

and Ad — o-i6A/p cos 0 , . . . (39) 

so that the rough estimates made here give a semi- quanti- 
tative explanation of the broadening and displacement of 
the lines. The apparent particle size can be estimated 
from both the apparent particle size and the displacement. 

4. Line broadening due to three-dimensional mistakes 
in AuCu 3 has been extensively investigated by Jones and 
Sykes (Sykes and Jones, 1936 ; Jones and Sykes, 1938), 
who have measured the apparent particle size for various 
reflections. These are given for a typical specimen in 
table II, most of the entries being average values from four 
or five photographs. It will be seen that e varies markedly 
with hy k and /, the extreme values, those for 100 and no, 
differing by a factor of two. The reason for the broaden- 
ing is fairly clearly the fact that the type of ordering is not 
the same throughout each single crystal, but ‘ changes 
step ’ occasionally, so that the probability of two cells 
separated by a distance t in the h k I direction being the 
same decreases with increasing t. The boundary between 
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regions of different types of order must be plane, as it 
would otherwise be unstable (Bragg, 1940), but there is no 
other general restriction on the mode of occurrence of 
mistakes. Wilson (1943) f“^2.s considered five different 

possibilities : 

(i) mistakes occurring on any plane, whatever its 
inclination to the crystal axes, 

(ii) mistakes occurring only on planes perpendicular to 
the edges of the unit cube «ioo> planes), 

(iii) mistakes occurring only on planes perpendicular to 

a face diagonal of the unit cube «iio> planes), 

(iv) mistakes occurring only on planes perpendicular to 

the body diagonal of the unit cube (<11 1> planes), 
and 

(v) mistakes occurring in such a way that gold atoms do 

not come into contact. 

The fifth of these resembles the second in some ways, as 
a little consideration shows that it is only on planes per- 
pendicular to the cube edges that mistakes can occur with- 

TABLE II 

Observed and calculated apparent particle sizes in a specimen of 

AuCu^ 


Line 

^ohs. 

^calc* 

% error 

100 

300A 

303A 

+ 1*5 

X 10 

100 

X09 

-+* 9*0 

310 

X 30 

131 

-h 9*3 

31 X 

X 43 

135 

— S '3 

331 

II5 

1 X 3 

— 3*6 

3x0 

123 

II4 

— 6*6 

330 

150 

152 

-f- X -3 


124 

123 

— 1-6 

410 -1- 333 

128 

137 

— 0-8 

4.11 

x6x 

157 

- a *5 

431 

IIS 

II 3 

— 3*6 


out bringing gold atoms into contact, but there is an added 
restriction. If, for example, the order is of type A, and 
a mistake occurs on a plane perpendicular to the ci axis, 
on the second hypothesis the order could change to B or 
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C or Z), while on the fifth it could change only to L>, The 
complete list of mistakes possible without bringing gold 
atoms into contact is 

on planes perpendicular to the a axis, A ^ D, B ^ C, 

on planes perpendicular to the b axis, B ^ D, A ^ and 

on planes perpendicular to the c axis, C ^ Z), A ^ B, 

Values of J{t) for each of these modes of occurrence of 
mistakes can be worked out by methods similar to those 
already employed for layer structures, but slightly more 
complex, because of the three-dimensional nature of the 
mistakes. If /?, k and I are all even or all odd 

J{t) = (Au + sCu)^ 

where Au and Cu are the scattering factors of the gold and 
copper atoms respectively- These are the reflections 
common to the ordered and the disordered state, and are 
sharp, in agreement with experiment, since J{t) is constant. 
If h, k and I are not all even or all odd the five possibilities 
of occurrence of mistakes lead respectively to the following 
values of J{fy and e : 

Variety of Apparent 

Mistake J{t) Particle Size 

(i) (Au — Cu)®exp{ — 4alf!/3<2} ' 

(ii) (Au — Cu)®exp{ — 4a/zlt|/3A/JV<3 } 3-v/iVa/2a^ 

(iii) (Au“Cu)^exp{ — 4a(A + A)|fl/3'\/ (zJsO^J } 3 •\/(2JV)<2/2a(/7 + ft) 

(iv) (Au — Cu)'‘*exp{ — 4a(ft+ft+i)!tI/3'\/(3N')<3 } 3V(3N’)a/2 a(ft+ft+Z) 

(v) (Au — Cu)“exp{ — aa(p + JSfa} JS[a/<x{p + q), 

where A, k, I are the positive values of the indices arranged 
in the order h > k ^ I, N — ->r P‘, a is the chance 

of a mistake occurring in a distance equal to a (the edge 
of the unit cell), and p and q are the indices of the same 
parity (that is, the two that are even, or the two that are 
odd, as the case may be). It is not possible to choose a 
so that any of these expressions gives quantitative agree- 
ment with the experimental values reproduced in table II, 
but the fifth comes close to it. The chief discrepancy is 
the 100 reflection, for which the theoretical value of e is 
a/oLip + o) = 00, so that the reflection should be sharp. 
It is actually the sharpest of the broadened reflections, 
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which suggests that a combination of two varieties of 
mistake, the fifth which gives good qualitative agreement, 
and the first to give a general fuzziness independent of 
indices, would give quantitative agreement. For this 
combination it is found that 

(vi) J{t) = (Au — Cu )2 exp { - \zoL'{p + q)/V^ 

+ 4«"/3]( l*IA) }, . . . ( 40 ) 

where a' and a" are the probabilities of the two kinds of 
mistake. This equation gives 


I 

+ g) ^ 2^ 
'\/Na 3<2 


(41) 


With this expression for e it is possible to get fair 
quantitative agreement with the experimental values. 
Apparent particle sizes calculated with cf! /a — 0*00303 
and /'^a = 0*00492 are given in column three of 
table II, and it will be seen that the greatest difference 
between the observed and calculated values is less than 


10 per cent. This agreement is very satisfactory, the 
difficulties of measuring line breadths being considered, 
but it does not seem impossible that some other mode of 
occurrence of mistakes would’ give even better agreement. 
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CHAPTER VII 


THE USE OF COMPLEX QUANTITIES 

I. In the earlier part of this book an attempt has been 
made to keep the necessary mathematics as simple as pos- 
sible, but a stage has now been reached where it seems 
desirable to introduce complex quantities. It is possible 
to proceed a little further without making explicit use of 
them, but already in the discussion of crystals with mis- 
takes it has been necessary to postpone one or two topics 
because elementary methods became very clumsy, and 
clumsiness is evident from the first in any attempt to dis- 
cuss diffraction by distorted crystals. It is, of course, 
impossible to give an adequate discussion of the nature 
and properties of complex quantities in the space available, 
but it is convenient to summarize the properties that will 
be used in the following chapters, and incidentally to clear 
up the difficulties postponed from chapters V and VI. 

A complex quantity is of the form 

z ^ X iy, . . . (i) 

where f ^ i and x and y are ordinary numbers or vari- 
ables. It is usual to call x the real part of z and y the 
imaginary part. A complex quantity may be represented 
in a plane by a line drawn froba the origin of length equal 
to r, the positive value of \^(x^ + and making an angle 
0 with the horizontal axis, where r cos Q — x^ r sin 0 = y. 
Since 

exp (z 0 ) = I + £0 — 0 ^ 2 ! —£ 0 ^/ 3 ! + 0 ^ 4 ^ + £ 0 V 5 * “ ”* 

= I — + 0V4- “ “ 0 V 3 - ^Vs* ” ***) 

= cos 0 + £ sin 0 , .... . ( 2 ) 

a complex quantity may be written 

z ^ X iy 
— r cos 0 + ir sin 0 
= r exp (£ 0 ). - 

85 


(3) 
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The sum of two complex quantities and -^2 is 

<2^1 + -s-g ^ + iy-^ + ^2 + ^V2 = + ^2 + Kyi + • (4) 

and their product is 

= (^1 + zjyi)(^2 + iy^) 

“ (^1^2 -yiy2) + ^‘(^ 13^2 v+ ^23^1). . (5) 

The product may be written more neatly in terms of r 
and 6 , for 

'S'1^2 = exp (z( 9 i) . ^2 exp (262) 

= exp {i{ 6 ^ + 62)}- • . (6) 

If the sign of i is changed wherever it occurs in a complex 
expression, the result is called the conjugate of the original 
expression. It is usually represented by the same letter 
with an asterisk affixed, so that if ;2' = + zjy, 

= X - iy. . . . (y) 

A quantity multiplied by its conjugate is purely real, and 
equal to the square of its absolute magnitude, for 

zz^ -= {x + iy)(x - ly) = 4 . ^,2 ^ ^2 ^ (gj 

The relations between a;, y, r, 0 and are shown 

diagrammatically in 
figure I for two values 
of z. In forming the 
conjugate of a com- 
plex expression it is 
not necessary to re- 
duce it to the form 
X + jy first. For in- 
stance, the conjugate 
of (A - iB)/(C + tJD) 
is (A + 2 B)/(C ~ iD). 
2. These properties 
of complex quantities make possible a compact expression 
for the amplitude scattered by a unit cell. This is repre- 
sented both in amplitude and phase by 

F = Efj exp {z7ti{hx^ + ky^ + Izj ) } 



Fig. I . — Representation of complex 
quantities. 
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where is the scattering factor and Zj are the co- 

ordinates of the 7*^ atom in the unit cell, and means that 

the expression is to be summed for all atoms in the unit 
cell. It is easily verified that the real part of this is equal 
to A and the imaginary part to JB (chapter V, paragraph i). 

2. In the calculation of the scattering factor from various 
imperfect structures it is necessary to calculate the mean 
value J of AA' + BB' and K of AB' - A'B, where A, B 
and A\ B' are the A and B parts of the structure ampli- 
tudes of two unit cells with a certain separation in the 
crystal, li F ^ A + iB, F' ^ A' -h iB\ 

FF'^ = + iB){A' - iB') 

= {AA' + BB') - i(AB' - A'B). . (lo) 

Thus J is the mean value of the real part of FF'* and — K 
is the mean value of the imaginary part. It is often con- 
venient to represent J — iK by the single letter Y. 

3. It is therefore important to calculate the mean value 
of expressions like exp (2,7thix), where u is fixed and x is 
some variable quantity, such as one of the 'co-ordinates of 
position within a crystal, or the displacement of a unit cell 
from its mean position. When x varies over a small range 
of values on either side of zero, an approximate mean value 
can be obtained by expanding the exponential in a power 
series : 

exp { 2 ,Ttiux) = I + zniux — ZTt'^u^x^ — /SfTtHii^x^/^ + •*'. 

If positive and negative values of x are equally frequent, 
the mean values of Xy and other odd powers of x will 
be zero, and 

<( exp {z7tiux)y — 1 — Z7chi^<(x^)> + Z7i'^u^<^x^y ~ ***, 

where the angle brackets indicate average values. Now 

exp (— ZTt^u^^x'^y) = I — Z7t'^u^<^x^y + Z7t‘^u^(j>c^y^ — ••*, 

so that, if the difference between <Cx^y/2 and can 

be ignored, 

< exp (zjziux)y => exp (— Z7t^u^<(x^y), 



88 


X-RAY OPTICS 


In certain cases this approximation is exact, but it is always 
fairly good for small values of x. 

4. If the probability of x having any particular value is 
given by/(^c), where the integral oi f{x) from — 00 to + 00 
is unity, the mean value of exp {zTtiux) is given by the 
integral 

r+<» 

<(exp (2,'jziux)'y = 1 /(^) oxp (zTiiux) dx, 

J —00 

which will of course depend on u. For the present, how- 
ever, let us forget about the interpretation of this expression 
as an average, and consider the function g{u) related to 
fix) by 

giu) =• ^fix) es-p izniux) dx, . . (n) 

/ 

where the limits are understood to be — 00 to +00 and 
fix) is any function of x for which the integral exists, and 
does not necessarily fulfil the condition that Sf(x) dx ^ 
The function g{ti) so defined is called the Fourier transform 
oi fix). It has the curious property that 

fix) = ^giu) exp (— ^Tzixu) du. . . (12) 

Except for a change in the sign of fix) is also the Fourier 
transform of giu). If both /and g are purely real they are 
Fourier transforms of one another, without this slight lack 
of symmetry in the relation. There have been numerous 
examples of somewhat- disguised Fourier transforms in 
earlier chapters ; for example, the amplitude of reflection 
from a small crystal (chapter III, paragraph 6) is the three- 
dimensional Fourier transform of a function with value 
unity inside the crystal and value zero outside, and the 
intensity of reflection from a layer structure with mistakes 
(chapter V, equation (29) ) is the transform of — iK^n- 

5. A rigorous proof of the relation expressed by equa- 
tion (12) is difficult, but depends on the following argument. 
From the definition 

giu) ==^fit)fx.-pi^ 7 ziut)dt, . . (13) 

since the substitution of t for x as the variable of integration 
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makes no difference. Let 

F(oc) = J^(w) exp (— 27zixii) du. . . . ( 14 ) 

Then 

Fipc) = iUf(jO ip>7tiut) dt\ exp ( — 2.7€ixv) du, 
or, on changing the order of integration, 

Fix) = /[J exp ~ x)u} du'\f(t) dt. . . ( 15 ) 

The inner integral, J exp {z7ti(t — x)u} du, is very peculiar. 
1£ t = X it has the value infinity, since exp (o) is unity, 
whereas if t:^x exp {2,7zi(t — x^u} continually oscillates, but 
has the average value zero. The inner integral is thus zero 
if t:^x, but 2it t — X it has an infinitely high and infinitely 
narrow peak. In the integration with respect to t, there- 
fore, the only region contributing anything will be an infi- 
nitely narrow strip containing the point t = x, and F{p^ 
will be equal to f{x) multiplied by the area under the peak. 
In various ways this can be shown to be unity, so that 

=/(^) • • • (16) 

The difficulties over which much trouble has to be taken 
are first, justifying the change in the order of integration, 
and secondly, showing that the area of the peak is actually 
unity. One method of approach is to consider 

F(x, X) = J^(w) exp (— zTtixu) du . ( 17 ) 

instead of F{x), and to show that F{x, X) tends to f(x) as 
X tends to zero. 

6 . Let us now return to the problem of the scattering 
of X-rays from an imperfect crystal. The rays scattered 
by the unit cell at the point rij^a, njy, n^c (the unit cell, 
for short) are advanced in phase with respect to those 
scattered by the cell at the origin by 

= ^7c{n^u 4- 7i^v + n^w) . . ( 18 ) 

where u, v, w are the co-ordinates of Q with respect to one 
of the points of the reciprocal lattice (chapter III, para- 
graphs I, 6 ). If is the value of the structure amplitude 
( 9 ) for the unit cell, the amplitude scattered by the 

7 
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whole crystal will be 

^ G(u, V, w) =27F^exp (zd^), . - ( 19 ) 

n 

where 27 means that the expression is to be summed for 

n 

all unit cells. The intensity scattered will be given by 
equation ( 19 ) multiplied by its complex conjugate, so that 

H(u, V, zo) = £I:F„F',, exp {i(d„ - 

n n* 

=■■ SZFnF^- exp {37 k[k(mi - n\) + ^{ 11 .^ - w') 

n n' 

+ zo(n^ - 713)]}. . (20) 

Note that in calculating H it is necessary to sum for all 
unit cells separately in G and G^ ; 

I = 27-F^F* exp {i{dn - <3J} == 27 |FJ 2 . ( 21 ) 

n n 

is something quite different from since it is independent 
of u, V, w» It is in fact the total intensity scattered by the 
crystal in the hkl direction, or the integral of H over one 
unit cell of reciprocal space. 

In order to obtain from (20) an equation suitable for the 
actual calculation of H, it is necessary to get rid of the 
separate values of F for each unit cell, which are not known 
individually, and replace them by some sort of average 
values, which can be calculated in ways similar to those 
used for layer structures. Clearly the average must be 
taken in such a way that the exponential factor in equa- 
tion (20) is not changed, or in other words, FJF'^r must be 
averaged for all values of n and n' that give the same values 
of 72 ^ — W2 “ ^3 ” ^3- equation (20) let 

n\ = 72 i + 772 ^, 7t^ = + ^2, n\ = 723 H- Tfl^, . (22) 

so that 

H{u, V, w) exp { — 2:772* {m^u + + m^w) } 

n m 

exp { ~\-m^ +m^w) }, (23) 

^ 771 

proper attention being paid to the limits of the summations. 



THE USE OF COMPLEX QUANTITIES 9I 

What does the sum in brackets mean ? Starting with any 
cell, designated by the numbers we are to find the 

cell m-^a away in the direction of the a axis, away in the 
direction of the h axis, and m^c away in the direction of 
the c axis, then multiply the structure amplitude of the 
first cell by the conjugate of the structure amplitude of 
the second, and add for all possible unit cells. For how 
many cells is it possible ? At first one might answer : for 
every cell in the crystal, but actually this is true only for 
TWj ^2 =^3=0. For = I, for example, a layer of 
cells all over one side of the crystal must be omitted, as 
there are no cells at a distance a away in the direction of 
the a axis. For = 10, a layer ten cells thick must be 
omitted. Consideration shows that in general the number 
of cells to be included in the summation is the number 
common to the crystal and its ‘ ghost ’ shifted a distance 
in the direction of the a axis, mj) in the direction of 
the h axis, and m^c in the direction of the c axis. If this 
volume is called the number of cells to be included 
is V^/U, where U is the volume of one unit cell. If, 
further, is the average value of the product 
for all unit cells related by the displacement m^ay m^by m^Cy 
then equation (23) becomes 

H{^y V, w) = exp {—2ni{mjU+ni.2V+mQw)}, (24) 

m 

This is the general expression for the intensity of reflection 
as a function of position in reciprocal space for a crystal 
whose unit cells are not all alike, and is the solution of the 
problem postponed in chapter V, paragraph 9. It is clear 
that is a generalization of the V{t) used in chapter IV, 
and of where and are the functions 

defined in chapter V. From the definitions it follows that 
changing the signs of the m"s leaves unchanged, but 
turns F,^ into its conjugate. A little care is necessary, 
therefore, in carrying out calculations, as the algebraic 
expressions for and F,„ do not always take account of 
these properties automatically. 
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The summation over m in equation (24) is to be done 
for all values of Wg, ^^3 which the product to be 
summed exists. Since the sign of the non-existence of the 
product is the vanishing of we may write the limits 
symbolically as ± 00 for each of the m’s, remembering to 
stop the summation when vanishes. If, w, v, w are 
not too big, and F,,, is a ‘ smooth ’ function of the w^s, 
the summations in equation (24) may be replaced by 
integrations, giving 

H{u, z;, w) = ( Uabc)-^ [ I 1 

J — QoJ — 00 J —00 

exp {— 27ii{ux/ci •+ "vy /h + wss/c)'^ dx dy dz, . (25) 

where x = m^a^y = /Wg^, ^ F(;», y, ^), Y(x,yy 

are F„j, with the appropriate substitutions. In a crystal 
containing mistakes Y{x, y, z) may decrease so rapidly that 
the variation of V{x, y, z^ can be ignored, and the intensity 
written 

Hiu, V, w) = (Uabc)-wi 5 iY(x,y, z) 

exp {- 27ii{ux/a + vy/h -1- wz/c)} dx dy dz. . (26) 

For a crystal without mistakes Y(x, y, z) is a constant, 
equal to and the intensity becomes 

H{u, V, w) = {Uabcy^F^nSV{^^y> 

exp {— 27zi(ux/a + vy/h + wz/c)} dx dy dz^ . (27) 

giving an expression alternative to equations (i5)-(i7) of 
chapter III. It is easy to verify that equation (27) gives 
the same results in particular cases as were obtained by 
elementary means. For instance, for a parallelepiped of 
edges X, F, Z, 

F(*, z) = {U/abc){X -\x\){,Y- \y\ )(Z - |^1 ), (28) 
and equation (27) gives on integration 
H(u, Vy w) 

^ ^2 sin^ (jiuX/a) sin^ (nvY /b) sin^ {tzwZ/c) 

(nuy ‘ (TtTjy * (tzw)^ ’ ‘ 

in agreement with chapter III, equation (7). 
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m 


7. Equation (26) may be used to discuss the intensity of 
reflection from a structure with a variable interlayer spacing, 
postponed from chapter VI. The values of and 
found there (end of paragraph i) give 

” Jm 

« exp ( — ZTtilg) 


H- 


— ~ a)^~^ exp (— 'zni.zlg) 

■<w * 


+ * • • + exp ( — 'znimlg)'\ 

— i^2[i — a + a exp (— . . • (30) 

by the binomial theorem. The m in this equation is what 
would be called in the notation of the present chapter ; 
since depends only on m^, and not on m-^ and 

H{Uy Vy w) drops rapidly to zero very rapidly with u and 
zv. If a is small i — a[i — exp ( — Z 7 cilgy\ is approximately 
equal to exp {— a[i — exp (— zjttlg)']}^ so that 

exp {- 7W3a[i - exp (~ ZTtilg)']}^ 
y(^) = exp {— a[i — exp (— zmlgjjz/c}, . (31) 

Equation (26) for the intensity of reflection as a function 
of Uy Vy w reduces to 

p+QO 

H{w) -= iy/V) Y{z) exp (- Z7tiwz/c)d{z/c)y (32) 

J — 00 

or, since Y{- z) ^ Y^{z) (page 91) and F/C/ - iV, 


H{w') 


= exp {\_—Z 7 tiw—cx.{x — exp {—Z 7 iilg)y\z/c}d{z/c) 

+ complex conjugate 

= exp (— + ZTziw}-'^ 

Hh complex conjugate] 

= iVjP^[{a(i — cos z:jzlg) + i{z7iw+0L sin Z 7 tlg)}~'^ 

+ complex conjugate] 
- cos 27 tlg ) 

— cos ZTclg)^ + (zTrry+a sin znlgY^' 
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Thus the maximum intensity of reflection is displaced by 
(jx/zn) sin znlg towards the origin of reciprocal space, and 
the integral breadth in a? is (a/a) (i — cos znlg). The dis- 
placement is in the direction to be expected, for the mistake 
increases the average c spacing, and hence decreases the 
distance of the maximum intensity from the origin of 
reciprocal space. Both the displacement and the broaden- 
ing depend on I, and for reflections such that Ig is integral 
they vanish, and become very small if Ig is nearly integral. 
The X-ray reflections will therefore be affected to very 
different extents, depending on the value of /. Hendricks 
and Teller (1942) have calculated in detail the ool reflec- 
tions of a partially dehydrated montmorillonite, which 
appears to consist of layers of thickness 10 A and 13 A 
with about equal frequency and a about one-half. They 
give no quantitative comparison with experiment, but state 
that the qualitative agreement is good. Equation (33) 
cannot be expected to hold for such a large value of a, 
but, since is of the form H{w) can be got by direct 
summation of equation (24) as a geometric series. 

8. The methods of chapter V, paragraph 9, were sufii- 
cient for expressing the profiles and integral breadths of 
lines on powder photographs of imperfect crystals. Never- 
theless, expressions in terms of complex quantities are some- 
times more convenient for calculation. If Y{t) = J{t) —iK(t) 
is the mean value of the product FF^ of the structure 
factors of two cells separated by a distance t in the hkl 
direction, it is readily found that the line profile is given by 

I{S) = U-^iV{t)Y{t) exp {znist) dt, . . (34) 

the apparent particle size is, if the maximum is not displaced, 

e = IJI - [KF(o)]-i;F'(i)y(i) dt, . (35) 

and the integral breadth is given by X/s cos 6 . For 
crystals in which Y(f) decreases rapidly with Jf, so that the 
variation of V(t) can be neglected, these reduce to 

I{S) = (y/U)^Y{t) exp {znist) dt . . (36) 

and £ = [F(o)]-i/F(f)^f. . . . (37) 
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9. Equation (26) may be used to discuss the distribution 
of intensity in reciprocal space for single crystals of AuCug 
containing mistakes. This was of course unnecessary in treat- 
ing line broadening in Jones and Sykes' powder photographs, 
but the actual distribution in reciprocal space can be investi- 
gated by the scattering from a single crystal (Guinier, 1945 ; 
MacGillavry and Strijk, 1946). Sufficient details have not 
yet been published to show whether or not the distribution 
of mistakes postulated to explain the powder photographs 
is correct, but it is of some interest to obtain expressions 
for the intensity of scattering for the various models. 

The values of Y(x, y, z) required follow readily from the 
values of J{t) for AuCus given on page 83. Since J(t) is 
the mean value of AA' + BB' for cells separated a distance 
t in the hkl direction, the real part of Y(x, 3;, z) will be J(t) 
with ht/^/N replaced by Xy kt/'\/N by y and It / by z. 
The imaginary part is zero, since K.{t\ is zero. Thus for 
the five varieties of mistake Y{Xy y, z) has the values 

(i) (Au - Cu)^exp {- 4a + z^)/3a}y 

(ii) (Au - Cu)^exp {-4a \x'\ /3<z}, 

(iii) (Au - Cu) 2 exp {- 4a( x'\ H- \y'\ )/2\/^ci}y 

(iv) (Au - Cu) 2 exp 4a( a;] + \y\ + |43r| )/ 3 a/ 3 ^}» 

(v) (Au — Cu)^ exp {— 2a( \x"\ + \y"\ )/a}. 

In expressions (ii) and (iii) x'y y\ z' represent Xy yy z 
arranged in the order |a:'| >\y'\ > and in expression (v) 
x" andy" are the co-ordinates corresponding to the indices 
of the same parity. The combination of mistakes of the 
first and fifth kinds gives 

(vi) (Au - Cu)^ exp {-[2a'( \x"\ + \y"\ ) 

+ 4 .(x"'\/(x^ + y^ + z^)/^]/^}, 

where oc' and a" are the respective probabilities of the two 
kinds of mistake. For a cubic crystal equation (26) becomes 

H(tiy Vy w) = Na-^ I 1 I Y{Xy y, z) 

J — 00 J “OO J — ^ 

exp {— Z7ti{ux + vy + wz)/a} dx dy dZy (38) 
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and the problem is now only that of performing the inte- 
grations with the values of Y(x, jy, z) given by expressions 
(i)— (vi). This is perhaps easiest for (iv) : 

H{u^ V, w) 

= iV(Au-Cu)2Q exp {-[4a 1^1 /Sa/s 

+ Z 7 tiux\/a}d{x / d) 

multiplied by similar integrals in y and z 



= iV’(Au— Cu)2j^j' exp {-[ 4 a/ 3'\/3 


+ complex conjugate 




multiplied by similar expressions in y and z 


= Ar(Au-Cu )2 


I. 


+ 


I 

^3 —Z 7 tiu_ 


L4a/3 VZ 4a/3 V3 

multiplied by similar expressions in y and 
= N{Am-CuY 

8a 8a 8a 


3 V 3 


3 V 3 


3 a/3 


(s-Ti)’"'”"’" (s-v;)**'””** . 

The integration for x" and y" in (v) is similar, giving 
H{u^ V, w) 


v (39) 


= iVCAu - Cu )2 


a 


a 


a^ + {Ttu'y a^ + {nd') 


(40) 


where u" and z;" are the variables in reciprocal space cor- 
responding to the indices of the same parity, and 
is the ‘ normal ’ variation of intensity of reflection with, 
the third variable ; it" is large only for w" o, and falls 
to inappreciable values for w" of the order of the lattice 
parameter divided by the particle dimensions. This is a 
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very peculiar distribution of intensity in reciprocal space ; 
the regions of high intensity are disc-like, with the plane 
of the disc parallel to the directions indicated by the indices 
of the same parity. Thus for the loo reflections the disc 
is perpendicular to the loo direction, for the oio reflection 
to the OIO direction, for the no reflection to the ooi direc- 
tion, for the 21 1 reflection to the loo direction, and so 
forth. The distribution of intensity in reciprocal space is 



Fig. 2 . — Reciprocal lattice for AuCuo. 


represented diagrammatically in figure 2 . The calculation 
for expression (vi) is more difficult, but it seems that the 
result would be qualitatively similar to that with (v). The 
regions of high intensity would be discs with the same 
orientations, but of appreciable thickness. Unfortunately 
the single photograph reproduced by Guinier (1945) is 
insufficient to show if the regions of high intensity are of 
this form or not. It does show, however, that the sharp 
reflections are surrounded by a region of diffuse scattering 
that cannot be explained by mistakes of the type treated 
here. It is probably an effect of lattice distortion, the 
copper atoms surrounding a gold atom being displaced 
somewhat from the lattice points. 

The integration for expression (i) is most easily done by 
transforming to polar co-ordinates with the axis passing 
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through the point w, t?, w. If ^ = \/{u^ + and 

r — ■\/{x^ 4- ^2 + ar^), equation (38) becomes 


roo rn r 

His) = iV 

J 0 J 0 J 


2 n 

exp {— [40C/3 + 2,7zis cos &\{r / a)} 
0 


= 2,7zN 


J 0 


2 sin 6 dd 


0 [40C/3 + cos0]^ 
8a 


{r/aY sin 6d(r/a) dd d^ 



As would be expected, the distribution is spherically sym- 
metrical, falling off ultimately as from its initial 

maximum of 277riV/8a^. The expressions for H corre- 
sponding to (ii) and (iii) can be obtained similarly, but 
are of considerable complexity, and the simpler cases 
illustrate the methods sufficiently well. An experimental 
investigation of the scattering in the neighbourhood of 
several points of the reciprocal lattice to show definitely 
if any of the models is correct would be of great interest. 
A preliminary investigation by Edmunds, Hinde and Lip- 
son (1947) and the work of Strijk and MacGillavry (1946) 
on reflections of the type hk^ are in qualitative agreement 
with model (vi). 
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CHAPTER VIII 


DISTORTED CRYSTALS 

I. In the crystals discussed so far the lattice has been 
perfect, or, in other words, the unit cells have been all of 
the same size and shape, though the contents of the cells 
may have varied throughout the crystal. Crystals in which 
the size and shape of the cell vary from place to place form 
an equally important but more difficult problem. The 
distortion of the lattice may be of several kinds. For 
example, at all temperatures above the absolute zero the 
atoms of a crystal are kept in continual thermal agitation. 
The agitation is best imagined as the superposition of 
a number of elastic waves traversing the crystal. These 
waves are similar to sound waves, but the amplitude is in 
general smaller and the frequency higher. The frequency 
is, however, low compared with that of the X-rays, and 
for purposes of calculation the atoms may be considered 
at rest, but displaced from their mean positions. This 
problem is treated in the next chapter. Other sorts of 
distortion are found in metals used for constructional pur- 
poses. Perfect crystals of most metals are soft and easily 
deformed, and metals in a useful form are generally in 
a state of internal strain. This may be produced mechani- 
cally in drawing or rolling, or ‘ chemically \ by treatment 
of an alloy in such a way that it begins to separate into two 
phases, but is not permitted to finish the process. Thus 
duralumin (about 96 per cent aluminium, 4 per cent 
copper) in equilibrium at room temperature would consist 
of a mixture of crystals of practically pure aluminium and 
of a compound AlgCu. As used, however, Cu or AlgCu 
separates in sheets within the aluminium crystals, but on 
lattices continuous with those of the parent crystals. The 
unit cells of the aluminium and the new structure will not 
fit together without a considerable amount of strain, and 
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the composite crystal is much harder than pure aluminium. 
In a sense the Cu or AI2CU sheets may be regarded as 
mistakes, but there is a variation of lattice parameter as 
well as of cell contents. Another example of a variation 
of lattice parameter within the individual crystals is given 
by alloys with compositions in the neighbourhood of 
Cu^FeNis- The alloys are in themselves of no particular 
importance, but it happens that the lattice parameter is 
approximately a sinusoidal function of position within the 
crystal, so that the theoretical treatment is comparatively 
straightforward and the interpretation of the photographs 
easy. 

2. The alloy Cu^FeNig has been investigated by Daniel 
and Lipson (1943, 1944). In equilibrium above about 
800° C. it consists of cubic crystals of one kind only, but 
below this temperature it consists of a mixture of cubic 
crystals of two kinds, differing in composition. If a speci- 
men in equilibrium above 800° C. is held at a lower tem- 
perature (say 650° C.) for a time insufficient to reach 
equilibrium, a state intermediate between that of perfect 
crystals of one kind and that of perfect crystals of two kinds 
is obtained, and can be maintained indefinitely by cooling 
the specimen to room temperature. In the intermediate 
state the lines on powder photographs are sharp, but are 
flanked by satellite lines, whose existence can be explained 
by the presence in the crystals of layers perpendicular to 
the crystal axes and alternately copper-rich and copper- 
poor. These will have a varying lattice parameter in the 
direction perpendicular to the layers, and the assumption 
that the variation is sinusoidal accounts for most of the 
details of the photographs. Since Fe, Ni and Cu have 
approximately the same atomic number the structure ampli- 
tudes of the cells will vary little, and may be taken as 
independent of position. 

Consider, therefore, a crystal approximately cubic, but 
with the lattice parameter varying sinusoidally in the 
X direction from a — d to a + d, the wave-length of the 



DISTORTED CRYSTALS 


lOI 


variation being L. Then the position of the unit cell 
will be given by 

+ {Ld/'ZTtd) cos ^ ^ 

Jn = 

since the distance in .the x direction between the origins of 
successive cells is 

+ i )<2 + {Ld/znd) cos {z7c{n^ + i)a/L'y\ 

— {n-jCL + {Ld/zna) cos (z 7 cn^a/Ly\ 

— a + {Ld / Z7ca)\cos, (znn-^a/L) cos {zna/L) 

— sin {zuzn^a/L) sin (zTta/L) ~ cos {zTtn-ya / 

which becomes 

a — d sin {z7tn^a/L), . . . (z) 

since a is small compared with L, and cos (zna/L) *= i, 
sin {zna/L) N zna/L. If A', h\ V (not necessarily integral) 
are co-ordinates in reciprocal space, the amplitude of 
reflection as a function of position in reciprocal space will 
be given by 

G{h\ k\ V) = {Z 7 ti{h'x,^ + k'y,, + Vz^ya} 

n 

= jF! 27 exp {z7zi[nj^h' ■+■ (Ldh' /zTca^) cos {z 7 tn^a/L) 

n 

+ nji'+nj.']} ..... ( 3 ) 

where F is the structure amplitude of one unit cell. As 
a function of k' and V this is normal, and has large values 
only for k' ^ /e, V = where k and I are integral. Its 
behaviour as a function of h\ however, requires special 
investigation. Let us consider therefore 

G{h') = jP^exp {zjci\nh' + {Ldh! /znaP) cos {z.Tzna / JO)^ 

n 

= FI7 exp {zTcinh'} .exp {zni{L6h' /zna^) eos (z7zna/L)} ( 4 ) 

n 

If d is so small that Ldh' /a^ is small compared with unity, 
the second exponential in equation ( 4 ) may be replaced by 
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p 

the first two terms of its power series : 

exp {zTtiiJLdh' / 2 , 710 ^) cos {znna/L)} 

= I + i{Ldh' /a^) cos {zTzna/L) 

= I i{Ldh* / 20^) {2Tdna/L') 

+ i{Ldh' / 20^) { — 27 tina/L}^ . (5) 

and equation (4) becomes 
G{h') — FU exp { 27 iinh' } 

4 - iFiLdh' /2a^)Zl C'xrp {zTtinih' + a/L)} 

n 

+ iF{Ldh' /2a'^)E exp {pjtinijh! — a/V)^. . (6) 

n 

The first of these sums has a maximum of amplitude pro- 
portional to F 2ith' = A, where h is any integer ; the second 
has a maximum of amplitude proportional to F{Ldh' /2a^) 
at A' a/L = hy where h is any integer ; and the third has 
a similar maximum at h' — a/L — hy where h is any integer. 
Thus, in addition to a strong reflection at the points of the 
reciprocal lattice, there are weaker ones at distances of 
iz a/L from it. The intensities are in the ratio 

\Ld{h - a/L)/ 2 a^Y‘ ^ i + a/L)/2a^'\^y 

or approximately {L 6 h/ 2 a^y^ * ^ * {L6h/2a^Y‘. The inten- 
sities of the extra maxima, or ‘ side-bands ’ as they have 
been called, increase as the square of the order of reflection. 

3. For comparison with this result it is interesting to 
work out the reflection from a structure in which the 
structure amplitude, instead of the lattice parameter, varies 
sinusoidally. 

If — nUy 

Fn - -^[1 + g cos {zTcna/L/ly 
where ^ is a fraction and L is the wave-length of the 
variation, the amplitude of reflection is 

G(A')=FZ'[i +g cos { 27 tna/L')'\ exp {2ninh'} 

71 

=F 2 J exp {2nmh'}+F(g/2)U exp { 27 iin(h' +a/L)} 

n n 

+F(g/2)Z! exp {27tm(h^ —a/L)}. (8) 

71 

As before, the three sums lead to maxima at h! — hy 
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h -iz a/ L, where h is any integer, but the intensities are in 
the ratio : i : so that the side-bands are the same 

in all orders. In particular, there are side-bands flanking 
the origin of the reciprocal lattice {h = 0 ), whereas these 
are of vanishingly small intensity for lattice parameter 
variation. Bragg and Lipson ( 1943 ) have devised gratings 
that illustrate the different diflfraction phenomena caused 
by mistakes, variation of structure amplitude, and variation 
of lattice parameter. 

4. Sinusoidal variation of lattice parameter in the x 
direction gives subsidiary reflections at =t a/L, In 
a cubic crystal there is nothing to distinguish the x direction 
from the y and z directions, and one would expect varia- 
tions in these directions also, either coexisting through the 
crystal, or taking place in separate regions. The investi- 
gations of Daniel and Lipson (1944) indicate that the latter 
is more likely. Points of the reciprocal lattice of the form 
h, ki I will thus be accompanied by three pairs of sub- 
sidiary reflections, the intensities being proportional to 
Those of the form A, k, o will have only two 
pairs, since one of the indices is zero, and those of the form 
hy Oy o will have only one pair. The separation of the 
subsidiary reflections from the main line in a powder 
photograph will not be a simple function of 6 , but will 
depend on the indices of the reflection. This may be 
seen as follows. (Compare chapter VI, paragraph z.) 
A reflection whose indices are h, ky I will reflect at an 
angle d given by 

{za/X)^ sin^ B ^ h'^ + + P. 

One whose indices are h a/L, k, I will reflect at the 
slightly greater angle 6 dS given by 

{za/Xy sin2 {6 + 66) ^ (h + a/Vf + + P-y 

or approximately 

(za/TCy sin B cos 6 66 = ha fLy 


66 


h X 

+hP + P)'zL 


.sec 6. 
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The separation of the side-bands is thus greatest for re 
tions of the form hoo, and least for those of the form 
If h, k and I are all non-zero and unequal there wi] 
three side-bands on each side of the main line, and ' 
overlapping with the main line and one another ir 
measurements of intensities difficult. 

5. Equation (8) for the amplitude of reflection : 
a structure with sinusoidal variation of structure ampli 
does not involve any approximations, but equation (6 
variation of lattice parameter depends on the assum] 
that Ldh! /a^ is small compared with unity. An 
treatment hardly comes within the scope of this book 
it may be remarked that the second exponential in equ 
(4) is represented exactly by the series 

exp {i(Lbh' /a^') cos {zTcna/L) } = J ^{Ldhf / a^') 

+ az /a^) cos {p.Tcna/U) 

— a J^{Ldh' /a^) cos {-ZTtn .za/L) 

— ziJJ^dh' /aP) cos (jzozn . / L) + 

where Bessel functions of the o^ 

order. Each reciprocal lattice point is the] 
flanked by a series of side-bands at distances of d= 
db za/L, =b 2^/Lj ••*, the intensities of the central refl< 
being proportional to J% and the intensities of the < 

being proportional to J% ••*. It is a property of 1 
functions that 

Jg+a/f +a7i + - = i, . 

so that the total intensity is the same as that from an i 
torted lattice. (There is only one central maximun 
two side-bands of intensity two of J% &c.) Eani^ 
Lipson (1944) find that the higher order side-ban 
appear, and that when all factors are taken into ac 
the intensities are in agreement with expectation. C 
the difficulties is ‘ extinction % a phenomenon that 
the intensities of the main lines relative to those of th< 
bands. The original paper must be consulted for full c 
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6. X-ray diffraction from aluminium crystals containing 
4 per cent by weight (1*75 per cent by atoms) of copper has 
been investigated by Preston (1938) and Guinier (1938). 
The phenomenon is somewhat complicated ; at room 
temperature layers of copper segregate on planes perpen- 
dicular to the axes of the parent crystal, but at higher 
temperatures (about 200° C.) the copper layei's * evapor- 
ate ’ and layers of AI2CU form. The AlgCu has not its 
usual structure, but one similar to that of fluorite (CaFg), 
apparently because this structure can fit into the aluminium 
lattice with less strain. The sheets of Cu atoms at room 

o o 

temperature appear to be about 45 A square, 200 A apart, 
and about two atomic layers thick. The layers of AlgCu 
produced at higher temperatures are considerably larger, 
and their diffraction effects sharper. No detailed treat- 
ment of the diffraction effects has been given in either case, 
but there is one very interesting phenomenon. The 
X-ray photographs show streamers from the ordinary 
reflections on one side only, in a way that indicates that 
the intensity of reflection in reciprocal space is not sym- 
metrical about the points of the reciprocal lattice. The 
h k I reflection has streaks running out towards 
+ I, k -h T, I \ and h, k, I + i \ but not towards 

h — T, k, I', h, k — I 'y or h, ky I 1. This is one of the 
few phenomena that are not effectively symmetrical about 
the points of the reciprocal lattice, and is at first sight 
rather puzzling. It is, however, qualitatively explained 
by a combination of varying lattice parameter and varying 
structure amplitude. Copper atoms are smaller than 
aluminium atoms, and contain more electrons, so that the 
scattering factor is high where the spacing is low, and vice 
versa. If both vary sinusoidally 

^ na {Ld/^Tcd) cos {2,7tna/JL)y 
and = jP[i + g sin {p!, 7 tna / Ly] ; 

in order to ‘ synchronize ’ the variations in the way postu- 
lated above the wave-length L must be the same in both 
8 
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and and one must vary as the cosine, the other as the 
sine. The amplitude of reflection as a function of h' 
is then 

GQi') = exp { 2 , 7 tih' x^/ a} . . (13) 

= FU\x + g sin (znna/Lyi exp {zuxih'\n + ijLd/'ZTtaF) 

n 

cos (zTtna/L)}}, 

which becomes, on expanding and neglecting products of 
small terms, 

G{h') = FU exp {zTtmh ' } 

+ tF[{Ldhy 2a^) — ^/ 2 ] 27 exp {z7tin{h' + a/L)} 

n 

+ iF[(Ldhy2a^)+g/z]2J {zjzin^h' — a/L)} . . (14) 

n 

The first of these sums gives a maximum of intensity 
proportional to F^ at h' =* /t, where h is any integer, the 
second gives a maximum of intensity proportional to 
F\Ldh' /za^ — g/z)^ for h' ^ h — a/L, and the third a 
maximum of intensity proportional to F\Ldh' /za^ + 
for ^ h -h a/L, It is clear that the inner maximum will 
be weaker than the outer, and will practically vanish if 
Lbh/a^ is approximately equal to g. This simple model 
gives merely a pair of extra spots accompanying the main 
reflection, but it is not to be expected that L will be constant 
throughout the crystal, and a proper treatment would 
undoubtedly give a continuous streak extending from the 
main reflection towards A + i, the intensity on the high- 
order side being greater than on the low-order side. 
Similar layers perpendicular to they and z axes will produce 
streaks running outwards towards k + 1. and / + i. 

7 * After suitable heat treatment (‘ annealing *) the spots 
or lines on X-ray photographs of metals are sharp, but if 
the metal is deformed by bending, twisting, rolling, draw- 
ing, hammering, &c., the spots spread out into Debye— 
Scherrer lines, and the lines become broadened. Clearly 
the original perfection of the crystal lattice has been lost, 
but there has been some controversy over the nature of the 
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imperfection. There have been three main suggestions : 

(i) that the metal is broken up into ‘ crystallites ’ so small 

(io“® to IQ”® cm. in linear dimensions) that dif- 
fraction broadening occurs ; 

(ii) that the metal is broken up into crystals (about io~'* 

cm. in linear dimensions) with differing lattice 
parameters ; and 

(iii) that the crystals of the metal remain fairly large (about 

cm. in linear dimensions), but are elastically 
distorted. 

The distinction between the second and third suggestions 
lies chiefly in the nature of the distortion. According to 
the second the distortion is merely a uniform expansion 
or contraction of the whole crystal, whereas the third 
includes non-uniform expansion or contraction, as well as 
twists and bends. 

A careful study of the variation of the line-broadening 
with wave-length, Bragg angle, and indices of reflection 
can decide between the suggestions. The first would yield 
a broadening increasing with 1 sec 6, and possibly varying 
somewhat with the indices of reflection, if the crystallites 
tended to have a more or less regular non-spherical shape. 
The second would yield a broadening increasing with 
tan dy but independent of Z and the indices of reflection. 
The third would yield a broadening increasing with tan 6 
and independent of X, but it would be independent of the 
indices of reflection only if the tensile strain were, on the 
average, independent of crystallographic direction. Most 
crystals (tungsten seems to be an exception) are not elas- 
tically isotropic, so that strain independent of direction 
would imply a variation of stress with direction. On 
general considerations of static equilibrium it seems likely 
that stress, rather than strain, would tend on the average 
to be independent of direction, with a consequent variation 
of the strain, and a dependence of the line-broadening on 
the indices of reflection. 
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On these qualitative considerations the experimental 
evidence is unquestionably in favour of the third sug- 
gestion. It is found that the broadening increases with 
tan d and not with sec b (Brindley, 1940 ; Stokes, Pascoe 
and Lipson, 1943 ; Smith and Stickley, 1943), is inde- 
pendent of the wave-length of the X-rays (the same refer- 
ences, and Lfipson and Stokes, 1943), and depends mark- 
edly on crystallographic direction (Wood, 4932 ; Stokes, 
Pascoe and Lipson ; Smith and Stickley). Though the 
experimental evidence makes it almost certain that lattice 
distortion is the main cause of the line-broadening, the 
theoretical problem of relating quantitatively the breadth 
of the lines and the properties of the crystal remains. It 
is as much a problem in the theory of elasticity as in X-ray 
diffraction, and possibly the nature of the glide processes 
by which excessive strains in the lattice are relieved must 
be taken into account also. No exact solution has been 
found, and perhaps each type of metal and each mode of 
deformation requires special treatment. It is, however, 
possible to give an approximate general theory. 

8- For calculating the line-broadening it is necessary 
to know Y(t)^ the mean value of FF* for two cells separated 
a distance t in the h h I direction (chapter VII, para- 
graph 8). In a distorted crystal the structure amplitude 
of a cell at {pc, y, z) will differ from that at (o, o, o) for two 
reasons : the cell is displaced bodily, and it suffers a rota- 
tion. It is hardly possible to evaluate y(i) without 
assuming a particular variety of distortion. It will, how- 
ever, be true in general that for large values of t the relative 
displacement of the cells will be large and random, and 
Y{t) will therefore vanish. For i small the displacement 
of one cell relative to another a distance t away will be 
ej^y the local value of the strain in the h k I direction, 
multiplied by t. The effect of the rotation on the phase 
of the scattered X-rays is small compared with that of the 
displacement, and may be neglected in a first approxima- 
tion. The structure amplitudes of the two cells will 
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therefore differ only through a phase factor 

exp { 2 , 7 tiSej,t}, . . ■ (15) 

where S ^ z sin d/X is the distance from the reciprocal 
lattice point in question to the origin.* The required 
function T(i) is then the mean value of 


FF* = F^ exp {zTziSef^t}, . . (16) 

If the fraction of the crystal for which lies between e and 
e -h de is y)h(e) de^ this mean value becomes 


Y(t) = exp {zTiiSet} de^ . (17) 

J —00 

and the apparent particle size is 


e = 


^-]-00 

[F(o)]-i Y{t)dt 

J — 00 

r r+" "I 

I 1 {^TiiSet} de dt 

J — 00 LJ —00 _ 

r+Qo n p+Qo 

I 1 *Aa(^) ^xp {z7zi{St)e} de 

J —00 LJ —00 


- s-^ 


d{St), 


= *S-V 7 .(q), (18) 

by the relations between Fourier transforms (chapter VII, 
paragraph 5). The integral breadth of the reflection is 

/? = yl/e cos 6 = z tan 6 /iph(o). . . (19) 

For a given reflection increases with tan 6 and is inde- 

pendent of A, as stated above, and 

^7i(o) = 2 tan d/p. . . . (20) 

Measurement of integral breadths thus gives the fraction 


* In a cubic crystal the components of the relative displacement 
along the crystal axes are eht.h/-\/Mf eh.t.h/’s/N^ eht.l/'s/Nt 
where = + The difference in phase between the 

rays scattered by the two cells is therefore 
Z'TriJk . e 7 ith / N + k.ehtk/‘\/N + l.ehtlf 's/JSl^/a 

— 2,TT 's/ Meht / a = 377(2 sin 9 /X)eht. 
The expression is, however, valid for non-cubic crystals (Stokes 
and Wilson, 1944). 
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of the crystal for which there is no strain in the hkl 
direction. In fact, it is fairly easy to show that when, 
proper allowance is made for other causes of line-broad- 
ening the line profile is — e cot B) within a proportion- 
ality factor. The negative sign merely indicates that 
places where the lattice is expanded reflect at lower angles 
than the ideal, and places where the lattice is contracted 
at higher angles. 

Equation (20) may be regarded as a complete, though 
approximate, solution of the diffraction problem. There 
remains the problem of evaluating ^hip) terms of the 
elastic, and perhaps the plastic, properties of the substance - 
One simple assumption is that ipjjp) is inversely propor- 
tional to the maximum strain in the hkl direction, and that 
this is given by the maximum stress, assumed independent 
of direction, divided by the value of Young^s modulus 
for the hkl direction. For a cubic substance this gives 


j 3 cot B == A + BH 


where 


4^max,hl> 

H = Ql ^ k ^ + ^ 2/2 / 2 ^ 2)/(^2 + ^2 + / 2)2 


(si) 


is the maximum stress and h2> *^44 s.re elastic 
constants. Stokes, Pascoe and Lipson (1943) have found 
that this equation gives quite good agreement with experi- 
ment, and that is approximately the same as the 

known values of the yield stress (Megaw et al., 1944, 1945), 
but Srmth and Stickley found only qualitative agreement. 
Stokes and Wilson (1944) have attempted to take shear 
stresses as well as direct stresses into account, using the 
approximation of paragraph 4, chapter VII, for the mean 
value of FF*. This gives 


Y{t) = <F 2 exp {zniSej,t}y = exp {- ZTi^SH^C^ef;^} 

and e - [l^(o)]-T Y{t) dt = S-\27t<,eiy)-i, 

/5 cot 0 = 
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Thus p cot d is proportional to the root mean square value 
of the strain in the hkl direction. On the assumption that 
the stress is on the average spherically symmetrical this 
leads to 

p cot d = ^/{A + BH) 

where A = 8 jr{(i 3 ii+ 2 s? 2 )<pfi> -45-12(512 +2yu)Of2>}\ 

B = - 325t{(5ii - 5ia)2 - AsliKpl2> : i 

KPi2y mean square values of the direct 

and shear stresses and H has the same value as before. 
The experimental results so far do not appear to be 
sufficiently accurate to decide whether (21) or (22) is the 
better approximation. 
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CHAPTER IX 


EFFECTS OF THERMAL MOTION 

I. The atoms forming a crystal are in a state of continual 
agitation about their mean positions. The period of vibra- 
tion of an atom is very small compared with the time 
required for taking an X-ray photograph, so that the 
photograph records the time average of the X-ray scatter- 
ing. The period of vibration of an atom is, on the other 
hand, large compared with the period of vibration of X-rays, 
so that at any instant the scattering is the same as if the 
atoms were at rest, but displaced from their mean posi- 
tions. For the interpretation of X-ray photographs it is 
therefore necessary to find the time average of the scatter- 
ing from a crystal in which the atomic displacements take 
on all values appropriate to its temperature. The calcu- 
lation falls naturally into two stages : (i) the determination 
of what functions of the atomic displacements are involved 
in the average scattering, and (ii) the determination of the 
average values of these functions in terms of the absolute 
temperature and the elastic constants of the crystal. 

2. Suppose that because of the thermal agitation the 
origin of the unit cell is displaced from its mean position 

of X/j^ fljCty 

=n^b +Y„,> ■ - • (0 

The effective structure factor of the cell will therefore be 
changed from F to 

~ F’ exp {z 7 ti{hXj^/a + kY^/b + . (2) 

It is here assumed that the unit cell is displaced as a whole, 
which is not strictly true for crystals with more than one 
atom per unit cell, but it will be seen later that the assump- 
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tion is justifiable. The function determining the scatter- 
ing from any imperfect crystal is the mean value over 
position in the crystal of for two cells separated 

by the translation mj?, m^c (chapter VII, equation 

(24) ). Thus the problem is to find the mean value of 

exp {z 7 ci\h{X^ ~ 

. (3) 

The average value of this over all cells in a crystal at any 
instant is equal to its time average for a given pair of 
cells with the required separation. It was shown in chap- 
ter VII that the mean value of exp (z-jc) is approximately 
oxp (— so that, if quantities belonging to the first 

cell are unprimed, and those belonging to the second cell 
are primed, the average value of equation (3) is 

= F^ exp {- Q, 7 i^C^h{X - X')/a + k{Y - Y')/b 

+ /(Z-Z')A]^>}, . (4) 

where the angle brackets represent a time average. As a 
matter of fact the result of using this approximation is here 
correct, though the proof of its correctness is difficult 
(Born, 1942—3). The expression to be averaged is easily 
rearranged to 

(hX/a +k Y/b +lZ/cY+{hX'/a Y'/b +IZ' /c) ^ 

-zQiX/a +k Y/b +lZ/c){hX' /a +k Y' /b ^IZ' /c). 

The first two terms are always positive, being perfect 
squares, and will have the same mean value. Let 

M = Z7t\(hX/a +kY/b + IZ/cyy,^ . (6) 

where M is necessarily positive. At first one might think 
that the average value of the third term would be zero, 
except for the special case = o, when it would 

cancel the first two. Cells close together in the crystal, 
however, will on the whole tend to move in the same way, 
the tendency diminishing with increasing separation of the 
cells. The average value of the third term will not vanish, 
therefore, but will start at — zM for the nCs all zero, and 




II 4 


X-RAY OPTICS 


decrease numerically to zero with increase in the nCs. Let 

Pm = ^K{hX/a -rkY/b + IZ/c) 

QiX' /a + kY'/b -t- lZ'/c)'y. . (7) 

Since this is symmetrical in the primed and unprimedl 
letters = P—m- Equation (4) becomes 

Y^ = exp (P^ - 2.M). . . (8) 

Since P ^ decreases from zikf to o as m increases, Y^ 
decreases from its initial value of P^, approaching the con- 
stant value exp ( — zTkf ) asymptotically. We have seen 
that a constant value of Y^ means a sharp reflection, and 
a value decreasing to zero means a broadened reflection. 
In the present case Y^ can be regarded as the sum of two 
parts ; a constant part P^ exp ( — zAf ), and a variable part, 
which it is convenient to write F^W^, decreasing from an 
initial value of P2[i — exp (— zAf)] to zero with increasing 
m. Equation (8) becomes 


(9) 


Y^ = P2 exp (- zM) + F^W^, 
where = [exp (P^) - i] exp ( - zAT), 

and the average intensity of reflection from the crystal as 
a function of position in reciprocal space is (chapter VII, 
equation (Z4) ) 

H{U^ Vy w) = Hl{Uy Vy W) + H ^{U y V y W) . (^O) 

where 


Pi(m, Vy w) = C/-ip2 exp ( - 2 ,M)ZV^ ^ 

m 

exp {— Z7ti(m^u + m<^v 4 - m^w)} 
and ^ ( I X ^ 

H^{Uy Vy w) = U~^F^ZV^W^ I ^ 

exp {— zjci{m^u + m^v + m^w)').j 

The first of these is the ordinary expression for the inten— 
sity of reflection from a perfect crystal, except that it is 
diminished by the factor exp (— zJVT). The second is like 
that for the intensity of reflection from a crystal containing 
mistakes. The scattering from a crystal in thermal agita- 
tion consists therefore of two parts, the first the same. 
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except for a reduction in intensity, .as that from a perfect 
crystal with stationary atoms, and the second a broad re- 
flection similar to that from a crystal with mistakes. Since 
is equal to the maximum of the broadened 

reflection will coincide in reciprocal space with that of the 
sharp. The total intensity, the integral of the sum of 
the two parts over one cell of the reciprocal lattice, is the 
same as that from an undisturbed crystal. In general 
decreases very rapidly with m, and the broadened part of 
the reflection may be written 

H^{u, V, w) - {V/Uabc)F^mW(x, y, z) 

exp {— zni{xu/a + yv/b + zzo/c)} dx dy dz . (12) 

where x = m-^a^ y = z = m<jC and the integrations ex- 
tend from - CO to :+ 00. Within a factor, therefore, 
is the Fourier transform of W. 

3. Logically the next step is to evaluate M. and in 
terms of the absolute temperature and elastic constants of 
the crystal, but a digression on their physical significance 
and qualitative behaviour is perhaps desirable. The inter- 
pretation of M. is the simpler, being times the mean 

square displacement of the cells in the hkl direction. This 
is easily seen for a cubic crystal, for then 

M = zn\{hX/a +kY/b + IZ/c^y 

/72 . . 72 

= + k'Y + rzyy 

= zn^S\{h'X -rk'Y I'zyy, . . (13) 

where h\ k', V, the indices of reflection divided by 

+ k'^ + 1% 

are the direction cosines of the hkl direction. The expres- 
sion in parentheses is the component of the displacement 
X, y, Z in the direction determined by h\ k\ or, in 
other words, the displacement in the hkl direction. For 
crystals of other systems the geometry is more complicated, 
but the result is the same. 
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Since the intensity of the sharp component of the reflec- 
tion is reduced by the factor exp (— aAf), temperature 
motion reduces the intensity of all sharp reflections. The 
reduction is on the whole greater for the higher orders, 
that is, those with large values of S = z sin 0 /A, but it is a 
smooth function of S only if the crystal is cubic or elastically 
isotropic. The mean-square displacement will be greatest 
in directions in which the restoring forces acting on the 
atoms are least, and the intensity of reflection from planes 
perpendicular to these directions will be lower than that 
from other planes. The effect is very marked in substances 
like graphite, in which displacements perpendicular to the 
layers are much easier than displacements in the layers. 

In a similar way it can be shown that is 47t^S^ times 
the mean value of the product of the displacements in the 
hkl direction of two cells separated by the translation 

m^c. This will depend not only on the ease of dis- 
placement in the hkl direction, but also on the relation 
between this direction and the translation relating the cells. 
Thus and M will be large for the same points of the 
reciprocal lattice, but for a given point P^ will be largest 
when the translation m-^a, ‘fnjb, m^c lies in directions for 
which the crystal is not easily deformable, since the cor- 
relation between the motions of the cells will be greatest 
for these directions, and it^will be least when the trans- 
lation lies in directions for which the crystal is easily 
deformable. The same is tme for W^n, ; it will be large 
for the translation lying along relatively undeformable 
directions, and small for it lying along readily deformable 
directions. The intensity of the fuzzy part of the reflec- 
tion is related to by equation (12), whence it follows 
that the variation of v, w) with w, v, w is the converse 

of that of with Wg, wig. Thus there will be streaks 
of comparatively strong intensity in directions for which 
the crystal is easily deformed, and regions where the scat- 
tered intensity is weak corresponding to directions for 
which the crystal is not easily deformed. This behaviour 
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is similar to that of the scattering from small crystals, dis- 
cussed in chapter III. The region of scattering is extended 
in directions for which the dimension of the crystal is least 
(face perpendiculars), and is contracted in directions for 
which the dimension of the crystal is greatest (comers). 
Small thickness corresponds to easy deformability and 
large thickness to difEcult. 

4. The evaluation of M and bears many resemblances 

to the calculation of the heat capacity of crystals. The 
motion of the atoms can be analysed in terms of ‘ normal 
co-ordinates % each of which is equivalent to a wave of 
definite frequency traversing the crystal. There are three 
times as many normal co-ordinates as there are atoms in 
the crystal, the number of * degrees of freedom ’ being the 
same whether the motion is regarded as made up of waves, 
or of motions of individual atoms. If there are more than 
one atom in the true unit cell, the waves divide into two 
groups, called ‘ acoustic ’ and ‘ optical h The former are 
of low frequency, with wave-lengths large compared with 
the cell dimensions, and the latter are of high frequency, 
with wave-lengths less than the cell dimensions. The 
number of waves in the acoustic group is 3 AT, where N is 
the number of unit cells, the remaining waves being optical. 
The energy of the crystal is equal to the sum of the ener- 
gies possessed by the waves. It is known from quantum 
theory that the average energy of a vibration of frequency v 
at the absolute temperature T is 




hv 




exp (hv/kT) — I 
i kT for T yy hv/k. 


+ ^hv 


(14) 


In this equation h is Planck’s constant and k is Boltzmann’s 
constant. If the number of waves per unit volume with 
frequencies between v and v + dv is Qy dvy the total energy 
of a crystal of volume V is 


E = VlEyQy dV, 


• (15) 
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the range of integration covering all frequencies. The heat 
capacity per unit volume is then == V~~^(dE/dT\. The 
well-known Debye expressions for E and follow from 
the assumption that the crystal behaves as an elastically 
isotropic continuum for frequencies up to a certain maxi- 
mum fixed by the requirement that the total number 
of frequencies should be three times the number of atoms. 

It is clear that this assumption applied to the problem of 
thermal scattering will obscure most of the interesting 
effects, since the mean-square displacement will be inde- 
pendent of direction, so that the weakening of the sharp 
maximum will depend only on and not on h, k, /, and 
the cloud of diffuse scattering around the sharp maximum 
will be spherically symmetrical and show no streamers or 
other peculiar effects. A detailed treatment of M and 
would, however, require an extensive knowledge of lattice 
dynamics, and we must confine ourselves to the Debye 
evaluation of M and a first approximation to the diffuse 
scattering. In order to keep the equations as simple as 
possible we shall assume that the axes are orthogonal and 
that there is one atom per unit cell. 

5. A single wave of frequency v traversing the crystal 
will displace the unit cell in the hkl direction by an 
amount varying sinusoidally with the time. Actually there 
are always three independent waves of a given wave-length 
— corresponding to the one longitudinal and the two trans- 
verse waves in an isotropic medium — and in the evaluation 
of M we may assume that one of these is responsible for 
the displacement in the hkl direction, and the other two 
are responsible for the displacements at right angles to it. 
The displacement of the unit cell in the hkl direction 
is then of the form 

A^CO&ZTtVty . . . (16) 

where is the amplitude of the wave of frequency v. 
To calculate M we want the mean value of as a func- 
tion of the absolute temperature, and this is obtainable by 
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way of the energy of the vibrating crystal. Since the mean 
value of cos^ is one-half, = ^<^Aly. The kinetic 

energy of the cell is — 2 , 7 th^^/uAl sin^ ZTcvt^ where fj, is 

the mass of the cell, so that the mean value of the kinetic 
energy is The mean energy of the cell is twice 

this, the mean values of the potential and the kinetic 
energies being equal, so that the total energy of the crystal 
due to this wave is 


E„ = 2,N:rcHyiA^,y, 

so that <ff> = = E,,/j^ 7 t^NfAV^, . . (17) 

The mean-square displacement of a cell in the hkl direc- 
tion is given by equation (17) summed for all values of v 
that are effective, so that 

M = 2.7z^S\S^y 

= 2^:jt^SKif<i^yvQ„dv 


S^VCF^Q 

6N/^j 


dv. 


(18) 


The factor of one-third comes in because only one-third 
of the waves are effective in producing displacements in 
the hkl direction. For an isotropic medium, as assumed 
in the Debye theory, = iZ7cC~^v^, where C is an aver- 
age value of the speed of propagation of elastic waves. 
The derivation of this equation will be found in any good 
textbook on heat, for example, in chapter XX of Heat and 
Thennodynamics by J. K. Roberts. Thus on the assump- 
tions of the Debye theory 


ZTtSWC^'^ j j 

where p is the density of the crystal. Since 


127lC 


"1 ' 

J 0 


Vv^ dv = 3iV, 


(19) 


ra - = 4^^- 

^ 


SN 


V. 


3 

m‘ 


3 


(20) 
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Substituting the value of from equation (14) gives 


M = 




m 

0 \j 


hv 


+ ^hv 


dv 


exp (hv/kT) — I 
or, with the substitutions 0 = hv.^/k, x = Q/T^ y = hv/kT^ 

Sh^S^T 


M = 




where 


T 

^(x) = - - 

^Jo€ 


+ 


X 


(21) 


2juk0^ 

y dy 

exp ly) - I 4 

The integral occurring in ilf(x) is not expressible in finite 
terms, but a table from which it may be obtained is given 
by Debye (1914) and in the Internationale Tabellen, p. 574. 
A selection of values of ^(x) is given in table I. The 
quantity 0 is the same as the ‘ characteristic temperature * 
occurring in the Debye theory of heat capacity. For a: less 
than unity (temperatures greater than the characteristic 
temperature) ip{x) is nearly constant, and 

. . . (az) 

2lJik 0^ 


within three per cent. For greater accuracy in this range 
the expression 

ijs { x ) = I + 5^2/36 . . . (23) 

may be used. 

6. Equation (ii) for the diffuse intensity v, %d) is 

more or less rigorous, but it is difficult to proceed further 
with it without making approximations. At not too high 
temperatures exp (P„,) in equation (9) for may be 

replaced by i + P^, giving 

TV,, = exp (- 2M)P^. . . (24) 

One can then evaluate in terms of the amplitudes of 
the normal co-ordinates, much as was done for JVI, and 
obtain v, zv) from equation (12). The whole process 

really amounts to finding the Fourier transform of 
and then finding the Fourier transform of the result, giving 
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again. The intensity at m, w is then proportional 
to with V satisfying the relation 

v^/C^ = u^/a^ + + w^/c^. 

The same result can, however, be obtained by a more 
elementary argument based on the treatment of a sinu- 
soidally-distorted crystal in the previous chapter, and as 


TABLE I 

lvalues of the function ^ipc) 


X 

o 

1 

2 

3 

4 

5 

tPQx ) 

1*000 

1*028 

I* 107 

1*333 

1*388 

1*571 

X 

6 

7 

8 

9 

10 



1*771 

1*984 

2*205 

3*433 

2*664 



space does not permit including both methods the simpler 
will be given. 

From the point of view of the X-rays each elastic wave 
traversing the crystal is equivalent to a sinusoidal distortion 
of the lattice. It has been shown that a sinusoidal distor- 
tion in the x direction of amplitude Ld/'ZJca and wave- 
length L produces a pair of satellite spots accompanying 
each reciprocal-lattice point at actual distances of ± x/L 
from it and of intensity {Ldh/ 2 ,a^)^ smaller. This result 
has now to be generalized in two ways. If the wave normal 
is not parallel to the x axis, the satellite spots will be dis- 
placed from the point of the reciprocal lattice along a line 
parallel to the wave normal, not along a line parallel to 
the X axis. Secondly, it is the component of the displace- 
ment in the hkl direction that produces the scattering round 
the point hkl of the reciprocal lattice, so that the effective 
amplitude of the wave is its actual amplitude multiplied 
by the cosine of the angle between the hkl direction and 
the displacement direction. A wave of amplitude A and 
length L will therefore produce satellite spots of intensity 
proportional to (nJTAS)^, where JT is the cosine of the angle 
between the displacement direction and the hkl direction, 
and S = 2. sin 0/A is the distance of the point hkl from the 
9 
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origin of the reciprocal lattice, the spots lying at distances 
of ± T./L, measured along the wave normal, from the point 
hkl. This reduces to the earlier result for a wave of 
amplitude Lb/2,7ca with its wave normal along the oc axis. 
For a given value of L there are three waves of different 
frequencies, speeds, amplitudes, and directions of displace- 
ment. The intensity at a distance ^ in a given direction 
from the point hkl will thus be proportional to 

7z^s\riAi + riAi + . • (^5) 

where quantities belonging to the three waves are distin- 
guished by the subscripts i, 2, 3. In calculating M we 
could assume that one of the three waves produced a dis- 
placement in the hkl direction and the other two at right 
angles to this direction, but here we have not assumed 
isotropy of the crystal, and all three must be taken into 
account. The average values of the squares of the A s 
are given by equation (17), so that the intensity is pto- 
portional to 

rrf£, riE^ 

zNij\_ vl vl vl , 

rriE^ riE^ riE^ 

ZN/IS^I Cf Ci Cl . 

where the a»*s and C’s are the frequencies and speeds of the 
three vibrations of wave-length L. This expression gives 
the intensity at 5 as a fraction of the intensity at the point 
of the reciprocal lattice, which is N^F^, or N^F^ exp ( — zJML) 
when allowance is made for the temperature motion. The 
three waves travelling in the opposite direction will con- 
tribute an equal amount, so one would expect the intensity 
at u, V, zo to be given by 
V, w) 


(26) 


NF^ exp ( - zM) 

fj. ■ s® 






L Cl 


Cl 


C| J ’ 


(»7) 


where = u^/a^ + v^/b^ + for an orthogonal lattice. 
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The correctness of the proportionality factor can be veri- 
fied by showing that the total diffuse intensity is equal, to 
the first order in M, to the intensity lost by the central 
spot. On the Debye theory there is a maximum value of v, 
independent of direction in the crystal, and hence the cloud 
of diffuse scattering is spherically symmetrical, and stops 
abruptly at a limiting value of s. For the purpose of check- 
ing the proportionality factor we may ignore the differences 
between the E*s and C^s, and take the total diffuse intensity 
to be the integral of equation (27) over the interior of a 
sphere of reciprocal space of radius v^/C. The sum of 
the J^’s will be unity under these conditions, and the total 
diffuse intensity becomes 

du dv dw 

= abc^^H^iu, w)d(M/a)d{v /h)d(w /c) 

Cvm/O 

— U I ds 

II Jo 

or, since s - i/L ^ v/C and NU = V, 

Y r . .SK^. ^E,. dv 

Jo 

= F^exp(- 2M).2MN, . . (28) 

by equation (19). The reduction of the total intensity of 
the central spot is 

NF^ - NF^ exp ( ~ 2M), . . (29) 

which is the same as equation (28) to the first order in Af. 
Agreement to higher orders cannot be expected, because 
of the approximations made in deriving equation (27). 

7. Equation (27) shows that the intensity of the thermal 
scattering increases with increasing order of reflection (S^), 
decreases with increasing distance from the points of the 
reciprocal lattice and for temperatures greater that 0 
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varies like ^S^Texp (- KS^T\ where Kis independent of T 
and (to the Debye approximation) the indices of reflection. 

It would thus have a maximum value for S'^T — i/K. The 
approximations made, which are equivalent to the expan- 
sion of exp (P^) to obtain equation (24), would hardly be 
valid for such large values of S^T, but it seems likely that 
the predicted behaviour is qualitatively correct, and that 
the intensity of reflection in the regions near the points 
of the reciprocal lattice does ultimately decrease with in- 
creasing temperature or order of reflection (Born, 1942--3, 
p. 317). The factor in brackets, involving the speeds of 
propagation of the vibrations, is in general strongly depen- 
dent on the direction of 5 and on the point of the reciprocal 
lattice concerned. It has been worked out for cubic 
crystals by Jahn (1942), but the expressions are too com*- 
plicated to be quoted and discussed here. 
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